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Abstract 



The aim of this thesis is to investigate quantum entanglement and quantum non- 
locahty of bipartite finite-dimensional systems (bipartite qudits). Entanglement is 
one of the most fascinating non-classical features of quantum theory, and besides its 
impact on our view of the world, it can be exploited for applications such as quan- 
tum cryptography and quantum computing. This circumstance has led to a growing 
interest and profound investigations in this area. Although entanglement and non- 
locality are ordinarily regarded as one and the same, under close consideration this 
cannot be taken for granted. The reason for this is that entanglement is defined 
by the mathematical structure of a quantum state in a composite Hilbert space, 
whereas nonlocality signifies that the statistical behaviour of a system cannot be 
described by a local realistic theory. For the latter it is essential that the correlation 
probabilities of such theories obey so-called Bell inequalities, which are violated for 
certain quantum states. The main focus of this thesis is on the comparison of both 
properties with the objective of understanding their relation. In terms of the anal- 
ysis of entanglement, recent methods for the detection are presented and discussed. 
Because of the fact that the correlation probabilities in general depend on the mea- 
surement settings it is necessary to optimise these in order to reveal nonlocality. 
This problem is solved for a particular Bell inequahty (CGMLP) by means of a self- 
developed numerical search algorithm. These methods are then applied to density 
matrices of a subspace spanned by the projectors of maximally entangled two-qudit 
states. This set of states has not only interesting properties with respect to our 
investigations, but also serves to visualise and analyse the state space geometrically. 
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Kurzfassung 



Diese Diplomarbeit setzt sich mit Verschrankung und Nichtlokalitat in bipartiten 
endlich-dimensionalen Systemen (bipartite Qudits) auseinander. Die Verschrankung 
ist eines der faszinierendsten nichtklassischen Phanomene der Quantentheorie, und 
neben ihrer Bedeutung fiir unser Weltbild findet sie Anwendung in der Quanten- 
kryptographie und der Quanteninformatik. Diese Tatsache hat zu wachsendem 
Interesse und ausgiebiger Forschung auf diesem Gebict gefiihrt. Verschrankung 
und Nichtlokahtat werdcn fiir gewohnhch als ein und dasselbe angesehen. Jedoch 
ist dies unter genauer Betrachtung nicht als selbstverstandhch hinzunehmen, was 
daran hegt, dass Verschrankung durch die mathematische Struktur eines Zustands 
in einem zusammengcsctzten Hilbertraum dcfinicrt ist. Nichtlokahtat hingegen be- 
sagt, dass das statistische Verhalten eines Systems nicht durch eine lokal-realistische 
Theorie beschrieben werden kann. Fiir letzteres ist wesentlich, dass die Korrela- 
tionswahrscheinlichkeiten solcher Theorien sogenannte Bell-Ungleichungen erfiillen, 
welche jedoch durch bestimmte Quantenzustande verletzt werden. Diese Diplom- 
arbeit dient insbesondcrc dazu, beide Eigenschaften miteinander zu vcrglcichen. 
Fiir die Untersuchung der Verschrankung werden aktuellc Separabilitatskriterien 
vorgestellt und diskutiert. Aufgrund der Tatsache, dass die Korrelationswahrschein- 
lichkeiten im AUgemeinen von der Messsituation abhangen, ist es notwendig diese 
zu optimieren, um Nichtlokalitat nachzuweisen. Dieses Problem wird fiir eine be- 
stimmte Bell-Ungleichung (CGLMP) durch cincn sclbstcntwickelten numerischen 
Suchalgorithmus gelost. Die besprochenen Methoden werden dann auf Dichte- 
matrizen eines Unterraums, aufgespannt durch Projektoren von maximal verschrank- 
ten Zwei-Qudit-Zustanden, angewandt. Diese Menge von Zustanden hat nicht nur 
interessante Eigenschaften im Bezug auf unsere Untersuchungen, sondern dient auch 
dazu, den Zustandsraum zu visualisieren und geometrisch zu analysieren. 
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Introduction 



Modern physical theories often contradict human intuition and this especially ap- 
plies to quantum physics. One of the most sensational quantum phenomena was 
recognised by Einstein, Podolsky and Rosen in 1935. In their seminal work [T] they 
presented a physical situation in which quantum theory seems to violate the princi- 
ples of relativity. Besides the fact that in such a situation best possible knowledge 
of the whole does not include best possible knowledge of its parts, particles seem to 
somehow influence each other at the point of measurement, even if these are space- 
like separated. The theory here predicts a new type of correlations which cannot be 
understood on the basis of a local realistic description. These correlations have be- 
come known as EPR correlation^. At that time, their conclusion was that quantum 
theory must be incomplete. It was thought that this peculiarity can be eliminated 
through a more fundamental theory based on local hidden variables. Such a formu- 
lation has never been found. Furthermore, the existence of such a theory can be 
revised experimentally via so-called Bell inequality tests. Experimental evidence of 
EPR correlations in 1982 by Aspect et al. (see [2]) has led to a growing interest 
in this subject and a new research field named quantum information has emerged. 
Many useful applications such as quantum teleportation, quantum cryptography and 
algorithms for quantum computers have been proposed. Entanglement and nonlo- 
cality are now being accepted and investigated by various physicists. 

Irrespective of intensive research within the last decades, there still remain many 
open problems from a theoretical point of view and one of the aims of this thesis 
is to outline several fundamental ones. In this thesis we investigate bipartite qudit 
systems, i.e. two d-dimensional systems, in order to gain deeper insight into the 
subject. The motivation for studying such systems arises from the fact that those 
are the most simple extensions of the two-qubit system to higher dimensions (An- 
other simple extension can be realised via multipartite systems, where the number 
of qubits is increased). At this point we might anticipate that investigations become 
more and more unfeasible with increasing dimensions. Systems with continuous di- 
mensions and/or a great number of parties are almost impossible to analyse with 
current mathematical tools. 

The thesis is organised as follows: In Chapter 1, we first provide an introduction to 
the mathematical framework necessary to study entanglement, followed by a closer 
examination of how it can be detected, quantified and classified. In Chapter 2, the 
focus is on the existence of local realistic theories. We discuss Bell inequalities in 
details, particularly, the CGLMP inequality. As it will be seen, in order to reveal 
the nonlocality of a quantum state the measurement setup, which can be expressed 
in the form of a Bell operator B, has to be optimised. In general, this is a very 
demanding subject and we contribute to the solution of this problem by presenting 
a self-developed numerical search algorithm. We end this chapter with further re- 
marks on the comparison of entanglement and nonlocality. The purpose of Chapter 
3 is to expose the properties of quantum states in a geometric context by applying 
entanglement detection criteria and our algorithm. The analysis of the state space 
of bipartite qubits first leads to a tetrahedron and motivated by its attributes we 
construct its extension for higher dimensions: the magic simplex. 

^or simply quantum correlations 
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1 Entanglement 



1.1 Basics 

1.1.1 Definition 

In this section we introduce the mathematical definition of entanglement using the 
most common notation of quantum physics - the Dirac notation, which utilises 
bra and ket vectors. As we do not go beyond bipartite systems, we restrict all math- 
ematical definitions to these in order to simplify our considerations. Consequently, 
all state vectors are elements of a composite Hilbert space IH.^^ = Ti^ ® Ti.^, where 
Ti^ and Ti^ are the Hilbert spaces of the two subsystems A and B. All operators 
acting on Ti."^^ will be elements of the corresponding Hilbert-Schmidt space , 
that is the set of bounded linear operators. We will utilise the Hilbert-Schmidt inner 
product defined by {A\B)jjg = Tr^A'^B), which induces the Hilbert-Schmidt norm 
ll^llji^^ = ^yTr{A'^ A). Usually the operations on the subsystems H"^ and are 
attributed to the fictive persons Alice and Bob, respectively. In 1935, Schrodinger 
recognised that EPR correlations are related to states G that cannot be 

written as tensor products of state vectors |^E''^) E H.^ and \^^) G H.^ . Consider 
an arbitrary state G Ti^^ . By choosing a basis d^"^)} of and of 

Ti^ any state of Ti^^ can be written in the form 

\'^n = ll^^A^)®\3") ' (1-1) 

with Cij G C and normalisation — 1- "^^^^ states that can be written 

in the form 

1^^^) = 1^^)®!^^) (1.2) 

separable. There is a class of states that cannot be written that way and we denote 
them as non-separable or entangled. A separable state can contain only classical 
correlations while an entangled state can contain quantum/EPR correlations. The 
definition implies that a state (of a subsystem) in general cannot be described by a 
state vector. It follows that open quantum systems need to be described by density 
matrices. Those also enable us to take into account decoherence and imperfect 
state preparations in experiments. For the density matrix formalism however, a 
generalised definition of separability for mixed states is required. If we suppose that 
a product state, regardless of whether it is pure or mixed remains separable under 
local operations and classical communication (see §1.4.2p we can infer that a state 
is separable iff it can be written in the form 

P^^ = Y.P^pf ® pf . (1.3) 

i 

with Pi > and YliiVi = 1- 

1.1.2 The two-qubit system 

We start with the two-qubit system which has the minimal number of degrees of 
freedom essential for entanglement. The qubit stands for a system with a two di- 
mensional Hilbert space Ti = with an orthonormal basis denoted by {|0) , |1)} 
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and is therefore the quantum mechanical counterpart of a classical bit. Neverthe- 
less, it should emphasized that there are major differences between the classical bit 
and the quantum bit. While a classical bit can either have the value or 1 a qubit 
can in principle store an infinite amount of information because of the infinitely 
many superpositions of |0) and However, this information is unusable since we 
cannot distinguish between two non-orthogonal states with a single measurement. 
Hence, we can only work with superpositions during information processing (quan- 
tum computations), the outcome however should be in an eigenstate of the measured 
observable. It should be mentioned that qubit systems are more than just simple 
examples of low dimensional quantum systems. This is because of their relevance in 
quantum optics (polarisation of photons in horizontal \H) or vertical \ V) direction) 
and experiments with spin | particles (spin up ||) or spin down ||)). As we have 
a two dimensional Hilbert space we need an operator basis with four elements to 
declare an operator. In many cases it is useful to work with the Pauli operator basis 
{1, (Ti, (T2, (TsjH, that is an orthogonal basis according to the Hilbert-Schmidt inner 
product {cFk\o'i)HS = 25fc«, {'^Wi)hs = 0, (1|1)_h'5 = 2. Any operator can be written 
in the form 

3 

O = aol + ^aiai, (1.4) 

i=l 

with all flj G C In the case of a density matrix O = p we write 

p=^(6ol + 6-5) , (1.5) 

which is equivalent to (11. 4p and has advantages with respect to the following restric- 
tions 

&o = 1 , (1.6) 
beR^, (1.7) 
< 1 . (1.8) 

The first restriction follows from the condition tr{p) = 1. The second is a necessary 
condition for p to be hermitian. Since p also has to be positive semi-definite we get 
the third restriction by computing det{p) = |(1 — ||fo||^). {det{p) > is a necessary 
condition for non- negativity, in our case it is also sufficient because of tr{p) = 1 
it is not possible for p to have two negative eigenvalues). As shown from this, the 
quantum state of a single qubit can be fully described by a three dimensional real 
vector b that lies within a three dimensional sphere with radius 1. Vector b is called 
the Bloch vector and the sphere, Bloch sphere. If the vector lies on the sphere, the 
state is pure; if it lies inside, the state is mixed. This follows immediately from 
det{p) = for = 1, meaning one eigenvalue has to be 0. This implies that the 
other eigenvalue has to be 1, and consequently the state is pure. 
In the following, we discuss systems of two qubits. To describe such systems we need 
a four-dimensional Hilbert space TC^^ = ® and we can choose, for example, 
an orthonormal basis of type {|0^) (g> |0^) , |0^) (g> jl-^) , |l^) (g) \0^) , |l^) (g) 
If there is no likelihood of confusion we can get rid of the indices A and B and the 

= (J J) = lo) <i| + Ii) = ( • 'd) = lo) <i| + i Ii) <^3 = (J = |0> (0| - |1> (1| 
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tensor products {|00) , |01) , |10) , |11)}. Before we discuss the Pauli operator basis 
for bipartite qubit systems in detail, we introduce the most well-known entangled 
states, namely the Bell states 

|^±> = i=(|01)±|10)) , (1.9) 

|<|.±> = i=(|00)±|ll)) . (1.10) 

The reader may convince himself that these states are indeed entangled, to be more 
precise, these states are maximally entangled as will be seen in our subsequent 
discussion of entanglement measures in §1.31 One (perhaps unexpected) feature is 
that they are all equivalent in terms of local unitaries. The application of a unilateral 
Pauli matrix {ai ® 1,(T2 ® 1,(^3 (g) 1} onto a certain Bell state yields another Bell 
state (up to a non-relevant global phase) 

(1-11) 
(1.12) 

(1.13) 

(1.14) 

In the next section we systemize the above mentioned attribute to get the gener- 
alised Bell states in higher dimensional systems. It is certain that properties and 
applications of Bell states could be further discussed. However, we continue with 
the analysis of our bipartite qubit system by examining the operators acting on the 
Hilbert space. As indicated before we now introduce the Pauli operator basis for 
the four dimensional Hilbert space. Any operator can be written in the form 

3 3 3 

O = al (g) 1 + ^ ttiCTi (g)l + ^6il(g)(Ji-F ^ Cijai (g) aj , (1.15) 

i=l i=l i)i=l 

with a,ai,bi,Cij G C Due to the fact that Pauli matrices are hermitian, it is 
apparent that a,ai,bi,Cij G M for all hermitian operators (for example observables 
or density matrices). Once again we change the notation slightly 

n,m=l / 





n : 








51 : 


^ 




0-2 ^ 


51 : 


<l>±> ^ 






51 : 







and specify the constraints for density matrices 

1 

4 



c=], (1.17) 



I 



(1.18) 

+ 11^1'+ 5^ tL<3, (1.19) 

n,m=l 

Tr (I*) (^|p) > V|^)g7^^^. (1.20) 



11 '^ti ^nm 
3 



The first restriction has to be fulfilled because Tr{p) has to be 1. Hermicity of p 
requires restriction two. The condition Tr{p^) < 1 implies the third one, while the 
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fourth condition is the non-negativity condition. As a conclusion, we can say that 
for two qubits any density matrix p G "Hjj^ can be fully specified by two real vectors 
r and s and a real 3x3 matrix with elements tnm- The vectors r and s determine 
the local characteristics of p related to the systems A and B, which can be seen by 
computing the reduced density matrices 

PA^TTBP=\{l+f-&) , (1.21) 

PB^TrAP=]^{l + s-a) . (1.22) 

It can be seen that there is no dependence on the parameters which can be 
regarded as correlation parameters since they reflect correlations of classical or EPR 
type. For this reason, the matrix with components tnm is sometimes referred to as 
a correlation matrix. For a given density matrix p all parameters can be obtained 
by the use of the associated algebra 

n = rr(ai® 1-p) , (1.23) 
Si^Tr{l®ai- p) , (1.24) 
tnm = Tr {an ' p) . (1.25) 



1.1.3 Bipartite qudit systems 

The investigation of quantum systems with few dimensions with the aim of gaining 
insight into the fundamental properties of quantum theory has been one of the most 
seminal concepts of quantum information and has led to interesting observations 
and countless applications of entanglement. The two-qubit system has allowed us 
to study entanglement in the absence of mathematical complexity caused by high 
dimensionality. However, within recent years, the focus has been on entanglement 
in quantum systems with more degrees of freedom. We now go beyond the familiar 
two-qubit system and concentrate on the entanglement of bipartite systems with 
arbitrary dimensions. 

The four Bell states have been useful for many quantum algorithms and seminal 
experiments. There is a very insightful way to generalise those maximally entangled 
Bell states onto Hilbert spaces H^^ — C^i^C^ with any desired d > 2. The analogue 
of the 1$+) state in C'' (g) C'' is a state of the form 

with an arbitrary orthonormal basis of Ti.-^ and of Ti.^. We briefly 

recall that Bell states are equivalent in terms of local unitaries. Consider the appli- 
cation of any local unitary transformation Ua ®Ub onto |fio,o) 

^ d-l 

Ua ® Ub 1^0,0) = -7^J2^A |s^> ® Ub |s^> • (1.27) 

^ s=0 

Since the transformations js'"^) = Ua I^"^) and \ s'^) = Ub arc basis transforma- 
tions giving the orthonormal basis {l-s'"^)} of H'^ and d-s'^)} of H'^, the resulting 
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state is once again a Bell state 

d-l 



s'^) . (1.28) 



If there are d"^ — 1 particular local unitaries that produce additional d"^ — 1 mutually 
orthonormal Bell states starting from |f2o.o)) we end up with an orthonormal basis 
of d"^ Bell states. Local unitaries with such properties are the Weyl operators, 
defined by the action 

Wk,i\s) =w''^'-^^\{s-l)modd) , (1.29) 
w = e*'"/'^ , (1.30) 

with k, I E {0, ...,d — 1}. The complete transformation on TC^^ is given by W^^i ® 1, 
producing d"^ orthonormal generalised Bell states 

\^k,i) = iWk,i ® 1) |l^o,o) • (1.31) 
We give the explicit verification of orthonormality 

{nm,n\^k,l) = ^ E ~ ^1 ® - /)mod d) ® \s) 

r,s=0 
d-l 

= ^ E Ui''^'"'^""'^"""\(r - n)mod d\{s - /)mod d) {r\s) 

^ d-l 

= - ^ w^^'~'^""^'~" \((s - n)mod d\ {s - /)mod dl 

s=0 c 

Onl 



k{s—l)—m{s—l) 



'd^ " 
s=0 

d-l 

"d 



s=0 

= Snl^km ■ (1.32) 

Weyl operators obey the Weyl relations 

Wj^iWk,^ = w^'Wj+k,i+^ , (1.33) 
Wli = ^k! = ^''W.k,-i . (1.34) 

Initially, the Weyl operators were rather contrived for the quantization of classical 
kinematics instead of the construction of a basis of orthonormal Bell states for 
Hilbert spaces Ti.^^ = ® C'^. In §3.21 we discuss how this has to be understood 
and how this concept can help us understand the symmetries and equivalences of 
quantum states. Since we are now able to construct a basis of Bell states, we continue 
with seeking a practical operator basis for Ti^f. For qubits, the Pauli operator basis 
has led to a simple presentation of density matrices via Bloch vectors. For qudits, 
we once again expand operators in the form 

d^-i 

O = aol + Y, ^i'^i , (1-35) 

i=l 
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with the dxd matrices {1, Fi, Td2-.i} forming an orthogonal operator basis . 
If we impose the operators {Fj} to be traceless we can fix the parameter = ^ for 
density matrices because of the constraint Trp = 1, just as we did for qubits 



P = ^1 + XI "i^i . (1.36) 



It remains open to show which {Fj} are beneficial for computations and parameter- 
isations. There are various candidates and we want to discuss two of them, namely 
the generalised Gell-Mann matrix basis and the Weyl operator basis. Both 
coincide with the Pauli operator basis for dimension two. We start with the defini- 
tion of the generalised Gell-Mann (GGM) matrices. For every dimension d we have 
d"^ — 1 matrices divided into three groups: 

1. f^i^iii symmetric GGM matrices 

Af = \j) {k\ + \k){j\ 0<j<k<d-l (1.37) 

2. ^ifti) antisymmetric GGM matrices 

Af = \j) {k\+t\k){j\ 0<j<k<d-l (1.38) 

3. {d — 1) diagonal GGM matrices 



0=0 



:i.39) 



The definitions imply that they are all hermitian and traceless. For proof of or- 
thogonality please refer to [0]. Due to hermicity of the GGM matrices all expansion 
coefficients have to be reals G M for hermitian operators, furthermore for density 
matrices we can again assign a real d"^ — 1 dimensional Bloch vector 

p = ^l + bA-A, (1.40) 

with Bloch vector = {{bl''} , {6^*^} , {b^}) and A = ({Af } , {A^'^} , {A'}), with 
restrictions for j, k, I given in the definitions of the GGM matrices. The vector 
lies within the Bloch hypersphere which precisely means ||6a II < ^{d-l)/2d, 
originating from the constraint Tr^p^) < 1 for density matrices. While for dimen- 
sion two, all vectors within the sphere result in positive semi-definite operators p, 
whereas in higher dimensions this is not the case. This means that there are areas 
within the sphere that are restricted for density matrices because of resulting p < 0. 
Unfortunately, until now no general expression or parametrisation has been found to 
avoid the holes within the sphere. Hence, this criterion has to be checked separately. 
The alternative, namely the Weyl operator basis is given by the matrices intro- 
duced in (11.291) leading to 

d-l 

H^;,; = ^e^^^ls) ((s + /) mod ci| k,l e {0, ...,d - 1} . (1.41) 
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Once again, we have the unity 14^00 = 1 and (f — 1 additional matrices forming a 
basis of mutually orthogonal operators. The proof of orthogonality is very similar to 
(11.321) and can be found in [6J. Hence, in the Weyl operator basis the Bloch vector 
expression of any density matrix is 

p = ^l + bw-W, (1.42) 

with a, (P — 1 dimensional Bloch vector h\Y and a vector W containing all opera- 
tors Wki except Woo- It is clear that one has to take into account the arrangement 
of the vector components. The main discrepancy between by/ and the foregoing 
Bloch vectors is that the components of hw can be complex and they have to sat- 
isfy -I = e^''''bli (to be understood modulo d) for hermicity of p. This is 
implied when comparing the definition f ll.4ip with hermicity = p. As before, 
the constraint Tr^p^) < 1 enforces the vector hw to lie within a Bloch hypersphere 
II^vkII < — ^)/d. Equal to the vectors 6a, not all vectors hw within this sphere 
are permitted, since some lead to matrices with negative eigenvalues. An expression 
in terms of the components hki has not yet been established. 

Finally we want to extend our Bloch-type operator expansion (11.351) to bipartite 
qudit systems. We can do this by generalising the expression (ll.lSp to 

(P^i d^~i d^-i 

O = al (g) 1 + ^ aiVi (g) 1 + ^ 6jl O + ^ djVi ® Vj . (1.43) 

i=l 1=1 *ii=l 



1.2 Detection of entanglement 

In §l.l.ll we have introduced the definition of separability and entanglement for pure 
and mixed states. Even though the distinction is well defined, in practice it is difficult 
to either find a separable decomposition or to prove that such a decomposition does 
not exist. For this reason, it is preferable to find operational criteria. As it will be 
shown, all known methods are either unfeasible or not sufficient to solve the problem 
completely. 



1.2.1 Reduced density matrices of pure states 

For all bipartite pure states, a necessary and sufficient criterion for separability arises 
through the form of the reduced density matrices. It is always possible to change the 
indexing of expression (II. ip . using one index running from 1 to ■ (Ib (dimension 
of T-Ca times dimension of 7^^) instead of two indices running from 1 to and 1 to 
ds, respectively. 



AB 



dA-dg 



:i.44) 



i=l 



with elements jz"^) G Ti.^ and G Ti.^ of some orthonormal bases. For a given 
state 1^^^) one can minimise the number of nonzero coefficients a, by the use of 

proper orthonormal vectors { ^"^^j and { 
the numeration starting with nonzero coefficients 7^ for i G {1, .., r} followed by 



} (see [3]). If we change the order of 
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= for i G {r + 1, ..,dA ■ ds} we end up with the Schmidt decomposition 



AB 



:B 



:i.45) 



i=l 



with r < min {c?^, rf^} and {Aj} known as the Schmidt coefficients. By computing 
the partial traces we reahse that {Aj} are the eigenvalues of the reduced density 
matrices 



i=l 

PB = Tr^(|vl>^^>(vl>^^|)=X^A. 



:i.46) 
:i.47) 



i=l 



Obviously, the resulting density matrices are diagonal and have rank r. Conse- 
quently, r is called the Schmidt rank. By comparing (11.21) with fll.45p we infer 
Schmidt rank r = 1 for all separable states and r > 1 for all entangled states. Since 
the only density matrices with rank 1 are pure states (they have only one eigenvalue 
Ai = 1) we come to the conclusion that all separable states result in pure reduced 
density matrices, while entangled states give mixed ones. Note that for practical 
purposes it is not necessary to construct the Schmidt decomposition of |\E''^^) before 
tracing out a system because the rank will always be equal due to basis independence 
of the trace. Hence, we have found a powerful tool to distinguish pure separable 
states from pure entangled ones since it is trivial to examine the mixedness of a 
density matrix. To summarise, a necessary and sufficient condition for separability 
of pure states is given by 



l"^^^) is separable ^ Tr^p^/^) = 1 



:i.48) 



1.2.2 Detection via positive maps 



As we know, in general a quantum state can only be fully described by a density 
matrix. Here, for a given state it is more difficult to determine whether it is separable 
or entangled. We begin with a rather theoretical consideration of necessary and 
sufficient conditions which was presented by the Horodeckis in 1996 f7]. 
Consider an operator fl G Ti-ljg and a linear map A : ~* '^hs- ^^Y 

the map A is positive (P) if it maps any positive operator in Ti.]^^ into the set 
of positive operators in 'H'jjg, in terms of mathematics if Tr{Q - Pi) > implies 
Tr(A(fi) • P2) > for all projectors Pi G Ti-j^g and P2 G Ti-j^g. Now, consider the 
extended linear map [A ® 1„] : TCffg^Ain T^'hs^-^u- Here stands for the set 
of the complex matrices nxn and 1„ is the appertaining identity. The map A is said 
to be completely positive (CP) if its extension maps any positive operator S G 
'hOfjs ® M-n into the set of positive operators Ti-^s ® -^n, i-e. Tr(S ■ Pi) > implies 
Tr([A (g) In] (S) ■ P2) > for all projectors Pi G n\js ® Mn and P2 G n]js ® Mn 
and all n G N. It might be a bit surprising that positive maps are not necessarily 
completely positive. All possible physical transformations are CP maps because the 
transformation of a density matrix has to result in another valid density matrix. 
For the solution of our separability problem we recognise the following. Consider 
a separable state p"^^ = J^iPiPf ® Pf ^ positive map A^ (not certainly CP) 
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inducing [A^ ® 1b] . The mapping of p"^^ gives 

[Aa ® 1b] (j^P^pf ® Pf ) =J2p^ (Aa (pf )) ® pf . (1.49) 



According to our assumption A^ being a positive map Aa (pf) only has non-negative 
eigenvalues and the same obviously holds for pf . All in all, we have a sum of positive 
operators weighted with Pi > yielding a positive operator. Hence for positive 
maps Aa the inequality [A^ ® 1b] {p"^^) > is a necessary condition for p"^^ to be 
separable. Furthermore, it can be shown (see [7]) that for every inseparable state p 
there exists a positive map Ap so that [Ap ^ 1] (p) < 0. The existence of such a map 
Ap is an impressive matter from a theoretical point of view, but for a given state 
p there is no recipe for how to construct such a mapping. Nevertheless, there are 
several approved positive maps for the detection of entanglement. One of them is 
the reduction map 

Ared{p) = lTr{p)-p . (1.50) 

The positivity of the map can be proven by 

(^1 ITr(p) - p\^) = Tr{p) (^|^) - (^| p |^) > V |^) e{\^)\ (^|^) = 1} . 

1 

Tr{p) is the sum of all eigenvalues (all positive or zero for positive operators) and 
(^'l p l^f) can maximally yield the largest eigenvalue, thus the map is positive. 
Therefore the reduction criterion [8j signifies that a separable state has to fulfill 
the inequalities 

PA®l-pAB>0 and 1®PB-PAB>0, (1-51) 
with the reduced density matrices pA and ps- 

Another method based on positive maps is called positive partial transpose 
(PPT) criterion [[9j. For this criterion the well known transposition is used to 
determine entanglement. For a given operator A, the transposition is defined by 

A = J2a,,\i){j\ , (1.52) 

It is obvious that this map is positive 

A is positive (^| = ^a^i \i) (j| ^) > V|^) G {|^) : (^|^) = 1} . 

id 

(1.54) 

Choose the vector |\E'') = |\1'*) and compute 

(^'1 1^') = (^*| A^ \^*) = J2 b) (^1 > . (1.55) 
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due to (11 .54^ this is always true for all |\&'), thus transposition is a positive map. 
As a result, we have found the positive partial transpose criterion. For a given state 
pAB the inequalities 

Pab = E ^'3\P^B ■ 1^) (^1 ® |J) (^1 > (1.56) 

ijkl 

and pl% = Y,{^j\pAB\kl)-\i){k\®\l){j\>0 (1.57) 

ijkl 

are necessary conditions for separability. The PPT criterion is stronger than the 
reduction criterion in all cases (see [8]). It should be mentioned that in the case of 
a two-qubit system TC^^ = ® or a qubit-qutrit system H"^^ = ® it is 
also a sufficient criterion (see [7]). Since PPT is not a sufficient criterion for higher 
dimensional systems, there exist entangled states with positive partial transpose. 
As shown in §1.4l this feature leads to a phenomenon called bound entanglement. 




1.2.3 Linear contractions criteria 



Additional operational criteria can be constructed from contraction mappings |lUj . 
The concept behind those criteria also issues from the extension of maps, but yields 
criteria that are independent from the previously discussed ones. Consider an oper- 
ator Q G Ti.jjg and a linear map A : Tij^g 'H^is- The map A is a contraction iff 
it does not increase the Hilbert-Schmidt norm, i.e. ||A (fi)||j^_5. < holds for 

all fl. Once again we define the extension of A analogous to the previous section 
by [A (g) 1„] : Hjjg ® — > ® A^„. A map A is a complete contrac- 

tion iff it is Hilbert-Schmidt norm non- increasing for all possible extensions n G N, 
||[A(8) 1„] (S)||^_5, < for all S G Ti-j^g (S> M.n- Now, for a separable state 
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J2iPiPf ® pf ^ contraction Aa we have 



HS 



[Aa®1b] J2piPt®pf 




HS 



[1.58) 



We have found a necessary condition for separabihty: 



I [Aa ® Ib] (p^'^) I 



HS 



< 1 



:i.59) 



Thus, for the detection of entanglement, the relevant maps are the contractions 
that are not complete contractions. For instance, transposition is such a map. In 
general there is no common method to construct the whole set of non-complete 
contractions. One noteworthy criterion based on the idea of contraction surely is 
the so called Matrix realignment criterion, which can detect PPT entanglement 
in some cases. The basic principle is to use a map 71 : Ti^f '^hs ^^^^ ^ 
contraction for all separable states within the composite Hilbert space Tifj^. While 
the contraction [A a ® 1_b] only affects a subspace, the matrix realignment map TZ 
affects the whole space Hfj^ in the following way 



7^(p 



AB\ 



^{ij\P 

ijkl 



AB 



\kl)■\^) {j\^\k) {l\ 



:i.6o) 



As we can see realignment interchanges the basis vectors {\j)} of Ti.^ applied on the 
left-hand side with {\k)} of Ti.^ applied on the right-hand side. This map satisfies 



|7^(p 



AB\ 



HS 
B 



< 1 



:i.6i) 



for all separable states p^^ = YliPiPf ® Pi (see [I2])- All states that violate this 
inequality are necessarily entangled. Some entangled states that slip through the 
PPT criterion can be detected via the realignment criterion and examples can be 
found in [TT] and [12]. 



1.2.4 Entanglement witnesses 

The last method for the detection of entanglement we discuss are entanglement 
witnesses. They originate from the Hahn-Banach theorem in convex analysis stat- 
ing that a convex set in a vector space (in our case this is the Hilbert- Schmidt space 
Tifj^) can be fully described by hyperplanes. In other words, there always exists a 
hyperplane that separates the convex set from the complement. Based on the fact 
that a hyperplane can always be expressed by a normal vector W (in our case the 
vector is an element of Ti^f and therefore a matrix) and the convexity of the set of 
separable states, we are able to comprehend a theorem stated by the Horodeckis [7] 
and B.M. Terhal [13]: A density matrix pab ^ '^//f is entangled iff there exists a 
hermitian operator W G 7i^5 with the properties 

TriWpAB) < , (1.62) 
TriWaAB) > , (1.63) 
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for all separable density matrices ctab- A hermitian operator W accomplishing those 
requirements is termed an entanglement witness. As a consequence, W has to 
have at least one negative eigenvalue and due to linearity of the trace, nonnegative 
expectation values on the subset of product states 

(^^1 (g)(^^|Vr|^^) ® 1*^) > \/ \^^) ® \^^) e n^'^ . (1.64) 

We can rank witnesses W by comparing the sets of entangled states they detect, 
that is Dw = {p\Tr{Wp) < 0}. An entanglement witness Wi is finer than W2 
iff Dyy^ C Dw-^. It is optimal if there is no other entanglement witness which is 
finer. Hence, an optimal witness WopUmai defines a tangent plane to the set of 
separable states which means that there must be at least one separable aAs with 

Tr{WoptimalO-AB) = 0. 




Figure 1.2: Schematic illustration of an optimal and an non-optimal entanglement witness 



Further information for the optimisation of entanglement witnesses can be found in 
. It is evident that a single witness does not detect all entangled states. We only 
know that for a given entangled state p there must be an appertaining entanglement 
witness Wp, an algorithm for the construction however has not yet been established. 
For this reason entanglement witnesses cannot be seen as a satisfying solution to 
the separability problem. Furthermore, it is also not clear how many witnesses are 
necessary to describe the whole set of separable states. Note that if the set is not a 
polytope then an infinite number of witnesses is required. The geometry of the set 
depends on the dimension of the considered space or subspace of . 



1.3 Entanglement measures 

The quantification of entanglement is another open problem. This is beyond the 
separability problem, since the purpose is not only to ascertain if a state is entangled 
or not, but also to quantify the amount of entanglement within it. This amount 
should capture the essential features that we associate with entanglement. Before 
we mention some possible measures we state some requirements on the attributes 
of a suggestive measure. Take into account that we restrict ourselves to bipartite 
systems Ti.^^ = with arbitrary dimension d. 
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1.3.1 The postulates 

In the first place, an entanglement measure E{p) is a mapping of density matrices 
into the set of positive real numbers 

p E{p) e R+ . (1.65) 

Normalisation is not necessary but reasonable and we therefore set 

E{\nk^i) {nk^iD = log2id) (1.66) 

for the Bell states 1^2^ ;). This should also be the highest value reachable, because 
we anticipate them to be maximally entangled due to maximal mixedness of their 
reduced density matrices pa = Pb = which is unique for pure states 

< E{p) < log2{d) . (1.67) 

By definition, the outcome of E{p) must be zero for all separable states psep 

EiPsep) = . (1.68) 

The most important requirement is monotonicity under LOCC. Entanglement 
cannot be created by local operations and classical communication ( §1.4.21 ). hence 
the contained amount has to be non-increasing under such transformations 

EiALOccip))<E{p) . (1.69) 

This implies that in the special case of local unitaries the outcome of E{p) will be 
invariant [I5] 

E (Ua ® UbpU\ ®ui^=E (p) . (1.70) 

These postulates are the only ones necessarily required and universally accepted. 
Some people tend to put further restrictions on entanglement measures. Some of 
them seem very natural and others can bring mathematical simplicity. One of them 
is convexity 

Convexity seems to be plausible as we cannot increase entanglement by mixing states 
and mixtures of entangled states can result in separable states (see Asymptotic 
continuity is another optional restriction 

lim Wpn-C^nWus^^^ \E {pn) - E {(Jr,)\ ^0 . (1.72) 

n— >oo n— »oo 

The interesting thing is that along with the next limitation called additivity we 
can obtain a unique measure for pure states [I5]. The additivity is separated in 
three types. We first consider partial additivity 

E(p«") = nE{p) , (1.73) 

which is the weakest addivity criterion. It only reveals that the entanglement content 
grows linearly with the number of pairs. We may expect that it is reasonable to 
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assume full additivity (following formula with an equal sign). However, it turned 
out that this would exclude some appreciated measure candidates. The weakened 
version is the subadditivity 

E{p^a) < E{p) + E{a) . (1.74) 

A more detailed analysis of the postulates and thermodynamical analogies can be 
found in [16], [17] and [18]. As mentioned before, the additional constraints of partial 
additivity and asymptotic continuity, inevitably require the measure to coincide with 
the von Neumann entropy S^n of the reduced density matrix p^^^ = Ttb/a {p^^) 
for a pure state p^^ 

EMp^"") = S.N (p^/^) = -Tr [pHog, [p^)) = -Tr [p^log, (p^)) . (1.75) 
This is the uniqueness theorem of entanglement measures ^19j . 

1.3.2 Measures based on distance 

The most intuitive measures are the ones based on distance. The concept is to regard 
the distance V of a state p G Ti.^^ to the closest separable state a E S C 'Hfj^ as 
the contained amount of entanglement 

E{p) = MV{p,a) . (1.76) 

crgcS 

It is apparent that those measures guarantee E{p) = for separable states. There 
are several types of possible distance functions V. They do not have to fulfill all 
criteria for a metric in a mathematical sense, but should meet the requirements of 
an entanglement measure. Let us consider the relative entropy of entanglement 
which was introduced by Vedral et al. [20] 



Er{p) = inf Tr (p {log2p - log2cr)) . (1.77) 

The distance used here does not meet the requirements of a metric because it is 
not symmetric and does not satisfy the triangle inequality. This is acceptable since 
relative entropy satisfies all criteria of an entanglement measure including asymp- 
totic continuity and partial additivity (see [21] )• Another distance is induced by the 
Hilbert-Schmidt norm and was investigated in [22]. The resulting measure is called 
Hilbert-Schmidt measure or Hilbert-Schmidt entanglement 

EHsip) = innp-a\\%, . (1.78) 

This measure has a different scaling than most of the other candidates because no 
logarithm is taken. For this reason, it is not common to normalise this function in 
the foregoing way. Up to now it has not been proven that monotonicity under LOCC 
is accomplished, meaning it is not clear if it is a good measure of entanglement (see 

mm)- 



1.3.3 Convex roof measures 



The idea of convex roof measures is to use entanglement measures for pure states 
Epure and generalize them to mixed ones in the following way 

Eip) = inf ^p.Ep..e (1^.) (^.1) , (1.79) 
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where the infimum is taken over ah possible decompositions of p = YliPi l^j) 
with YliPi — 1 Pi — 0- The decomposition yielding the infimum is then said to 
be the optimal decomposition or optimal ensemble of p. It can be shown that 
if the utilised measure for pure states Epure is monotonous under LOCC then the 
induced convex roof measure has this property too (see [5]). The first measure of 
this kind was entanglement of formation Ep with the von Neumann entropy of 
the reduced density matrices as a measure for pure states Epure = E^]\f introduced 
in ([TTSD 

EM = ^ inf Tp^E.^ (v&.l) . (1.80) 

When this measure was introduced by Bennett et al. |24j , the notion was to quantify 
the amount of pure Bell states that is needed per copy to construct p. In other 
words, copies of p can be prepared out of a minimum of Ef{p) ■ N Bell states 
and (1 — Ef{p)) ■ N separable states. Ep{p) can then be regarded as the contained 
amount of entanglement. The naming "entanglement of formation" is referable to 
this concept. 



1.4 Bound entanglement 

In the previous sections we frequently used the term local operations and clas- 
sical communication (LOCC). Nevertheless, we neither explained its meaning 
nor given an adequate mathematical definition. This section will serve as an in- 
troduction to this issue. After introducing the class of quantum operations, we 
consider the restrictions for the class of LOCC. Consequences for the purification 
of entanglement will be briefiy discussed, enabling us to justify the term bound 
entanglement. 



1.4.1 Quantum operations 

We investigate the fundamental quantum operations (see [I9], [25]). Our objects of 
interest are the maps that transform a given state p G 'H\jg into another p' G 'H'\jg 

^■.n],s^ nls , (i-si) 
P = Hp) ■ (1-82) 

As far as we know, any A is a combination of four elementary linear map^: 

• Adding an uncorrelated ancilla a G Ti-^s original quantum system 
in the state p G H^g 

Ai(p) = p®cr. (1.84) 

• Tracing out part of the system in the state p G Ti\is ® '^hs 

A2 : nhs ® nls ^ nj,s , (l-85) 

A2(p) = Trap • (1.86) 
Here Tr2 is the partial trace over the Hilbert-Schmidt space Ti-^g. 



^It is true that the first two types of transformations are not "real physical" operations, but they are the 
conversion of experimental operating principles into mathematical diction. 
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Unitary transformations of a state p G Ti}, 



A3 : n]js® n 

Asip) = UpU^ , 



HS y 



(1.87) 
(1.88) 



with a unitary operator U G Ti-jj^. 
• Measurement of an observable A E H 



HS 



A4 : n 



HS 



n 



]lS ' 



(1.89) 



Up) 



(1.90) 



t 



) 



Here {Mj} is tlie set of the spectral projectors associated with the eigenvec- 
tors of the observable A. It is useful to classify two types of measurements: 
Non-selective measurements, where we work with all outcomes of the mea- 
surement and selective measurements, where we filter outcomes. Non- 
selective measurements then yield YliM^Mi = 1, and selective measurements 



The point of the matter is, that the four maps themselves and their compositions 
are all completely positive. According to Choi's theorem (see [26]) any map of this 
type can be expressed in the form 



with p G Ti-lfg (Hilbert-Schmidt space of dimension n x n), A(p) G Ti-jjg (Hilbert- 
Schmidt space of dimension mxm) and complex matrices {^} of dimension mxn. 
This expression is the Kraus representation of A and the matrices {Vi} are the 
famous Kraus operators which obey Ei ^i^i — 1- 




(1.91) 



(1.92) 
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1.4.2 Class of LOCC 



The concept of "local operations and classical communication" stems from 
quantum communication theory and appears in quantum teleportation, quantum 
cryptography and distillation protocols. Consider a source and two distant parties 
A and B, commonly called Alice and Bob. Under realistic circumstances they are 
able to perform arbitrary quantum operations acting on the particular local Hilbert 
space Ti"^ and and to communicate via classical information. They can neither 
exchange their quantum systems nor perform "global" transformations involving the 
entire composite Hilbert space Ti^^ = Ti.^ (E) Ti.^. 



quantum channel / / / \\ \ quantum channel 

n iU ) w 



. ,. classical communication channel _ , 

Alice ► Bob 



Figure 1.3: Schematic illustration of the LOCC situation 

The transcription of this concept into restrictions for the quantum operations §1.4.11 
implies that a LOCC operation is a map of the form 

A.bM^ (1.93) 
rr[Ei(-4i®Bi)p(A*®B,') 

with product Kraus operators {Ai ® Bi} where Ai acts on Alice's Hilbert space Ti.^ 
and Bi acts on Bob's Hilbert space TC^ . The product Kraus operators reflect that 
the quantum operations only act locally, while the bilateral dependence of {Ai ® Bi) 
on i reflect that both parties can arrange their actions (they can be classically 
correlated). While the form of fll.93p is comprehensible, the constraints for the 
operators {A^ (g) Bi} and their relations can be quite complex. Those depend on 
the considered communication class, which can be "no communication", "one-way 
communication" or "two-way communication" (see [5],[l9],[25]). The form of (11.931) 
induces the definition of separable density matrices (II. 3p . 



1.4.3 Distillation and bound entanglement 

Outstanding procedures like quantum cryptography or quantum teleportation re- 
quire pure Bell states. In practice however, we cannot completely neutralise the 
interaction with the environment that causes decoherence. In the case of the two 
distant parties Alice and Bob, the linking quantum channel is therefore said to be 
noisy. We expect that the occurring state p"^^ is no longer a pure Bell state, origi- 
nally emitted from the source, but a mixture. We characterise the state p^^ by its 
fidelity F 

F = Tr (|fifc,;)(fifc,^|p^^) , (1.94) 

which can be regarded as the remaining content of 1^2^;) {flk,i\ in p^^. In order to 
reconstruct an almost pure Bell state with any desired fidelity close to 1 via LOCC 
operations we have to execute a so-called distillation protocol. The first protocols 
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for qubits were introduced by Bennett et al. [27] . These are the recurrence, hashing 
and breeding protocols. The hashing and breeding protocols were later generalised 
for qudits by Vollbrecht et al. [28]. These will not be discussed in detail, but it 
should be pointed out that they all have one thing in common: They use several 
copies of p^^ and local filtering to obtain a smaller number of nearly maximally 
entangled pure states. For every protocol there is a lower bound for the fidelity 
Flb in order to work successfully. It is a challenging task to find more universal 
(ones with a lower F^b) and faster protocols. Irrespective of this intention, there 
will always be entangled states that cannot be distilled. The conventional term 
therefore is bound entanglement (see [29]). We have already mentioned that in 
higher dimensional systems there are entangled states which cannot be detected 
via the PPT criterion. Let us consider how such a state behaves under a LOCC 
transformation. According to the assumption p G Hfj^ being an entangled density 
matrix with positive partial transpose we claim 



(^|p|$) > v|^) G n 

(^|p^^ 1^) > V|^) G ?i 



AB 



AB 



A LOCC transformation of p then yields 



Tr 



E,(A,®5,)p(4®5l; 



:i.95) 

[1.96) 



(1.97) 



We know that transposition changes the sequence of operators according to (LMN)'^ - 
N'^M'^L^ and in our case hermitian conjugation is simply the combination of trans- 
position and complex conjugation = {L*)^ . Partial transposition on Ti^^ then 
results in 



p = iV 5^ (A, ® Bl)p {A\ ® {B*y) 



(1.98) 



wherein is a positive real number A^ 
us investigate positivity 



1/Tr j:^iA^B,)pi4®Bj 



> 0. Let 



(vl/| p |vl/) = AT^^ {A, ® B*)p^- (4 ® (i?*)t) 

i 

>0 V I *i>6-K -4^1X311 

>0 V|^)g7^^^. 



(1.99) 



This proves that p is a positive operator. The conclusion is that a LOCC trans- 
formation cannot transform a PPT density matrix into a density matrix that is 
non-positive after partial transposition (NPPT). We have established that a pure 
maximally entangled state (NPPT) cannot be distilled from a PPT state, even 
though it contains some entanglement. This justifies the term bound entangle- 
ment. The existence of NPPT bound entanglement is controversial and a subject 
of recent research. 
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2 Bell inequalities 



2.1 Local realism versus quantum mechanics 

Since the verification of (special) relativity principles most physicists have believed 
that any fundamental theory is in compliance with local realism. 
Realism is the assumption of the existence of definite values for all possible ob- 
servables. That is, at each point of time these values genuinely exist, whether we 
measure them or not. "Ideal" measurements with no or marginal disturbance could 
therefore reproduce these pre-existing values. As a direct consequence, realistic the- 
ories are non-contextual. 

Locality reflects the key consequence of relativity, that is all interactions between 
distant objects are limited to the speed of light. Space-like separated objects are 
therefore independent of one another. 
Local realism is the unification of locality and realism. 

As obvious as these assumptions may seem, in 1935 Einstein, Podolsky and Rosen 
showed that quantum mechanics rejects these principles (see [1]). At this period of 
time, their conclusion was that quantum mechanics must be incomplete, meaning 
there must be elements of reality that do not appear in the theory. An introduction 
of such elements however, should restore local realism. Due to their hiddenness 
they have been called local hidden variables, regardless of whether they are hid- 
den in principle or just in experiments at that time. During this period, evidence 
had not been found yet, demonstrating that all local hidden variable theories 
(LHVT) are incompatible with the predictions of quantum mechanics in particular 
experiments. In 1964, John Bell derived that the correlation expectation values of 
local hidden variable theories fulfill inequalities that can be violated by quantum 
mechanics (see [30]). Moreover, these Bell inequalities enable us to experimen- 
tally revise the validity of either quantum mechanics or a local realistic theory. 
Before going into further details, it is essential to briefly reconsider why and under 
which circumstances quantum mechanics is said to be non-local. Nonlocality can be 
ascribed to entanglement and the measurement problem. This can be best under- 
stood by considering the Bohm version of the EPR situation (see [31] )• Here, two 
spin- 1 particles, for example electrons or protons, interact and their spins are in a 
maximally entangled singlet state |\E'~) = -^(iTi) — lit)) afterwards. Then both 
particles fly off in different directions freely while the spins remain in the singlet 
state. Subsequently Gz spin measurements are performed on both particles. 



According to the measurement postulate of quantum mechanics the spin will no 
longer be a superposition of up |t) and down ||) though either one or the other. 
The reduction of the state vector to one of the eigenvectors of the observable Oz of 
one particle forces the other particle's spin Oz to be anti-correlated. Here appears a 
new type of simultaneousness. The problem here is that if one measures Oz on one 



position of interaction 




45(in)-UT)) 



Figure 2.1: Schematic illustration of the Bohm-EPR situation 
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side, the particle's spin on the other side is affected instantaneously, irrespective of 
how far they are apart. If the wave function is regarded as a real physical object this 
means that causality is violated, since one particle affects the other with superlumi- 
nal velocity. When we say the measurements A and B are performed simultaneously 
then according to relativity there are different inertial frames of reference, in which 
A happens before B and vice versa. Fact is, quantum mechanics predicts that the 
spins are always perfectly anticorrelated. This phenomenon can be interpreted in 
different ways. When we accept the wave function as a real existing physical object 
we have to reject locality, since the reduction of the wave function is a non-local pro- 
cess. In the conventional version of the Copenhagen interpretation where the wave 
function is not a physical object of reality but rather a mathematical tool, realism is 
rejected in order to preserve localit}0. In the Bohm interpretation locality is given 
up to maintain realism. What can be concluded from those interpretations is that, 
if nature really behaves quantum mechanically, then it is not local and realistic at 
the same time. 



2.2 A Bell inequality for two-qubit systems: CHSH 

As previously mentioned, local realistic theories are in contradiction to quantum 
mechanics in certain experiments. In order to show this, we derive a Bell in- 
equality which holds for any local realistic theory but can be violated by quantum 
mechanics. This will be the famous Clauser-Horne-Shimony-Holt inequality 
(CHSH), which is a modification of the original Bell inequality permitting experi- 
mental revisal. 



2.2.1 Derivation of the CHSH inequality 

We want to describe the Bohm-EPR situation with a local realistic theory. In order 
to do this, we introduce the parameter A, which represents a hidden variable or a set 
of those. The spin of a particle then depends on the measured direction represented 
by a vector a and the parameter A. We know that a spin measurement on a particle 
has only two possible outcomes. Without loss of generality we assign the values -|-1 
and —1 to them. In a local realistic theory the observables of distant parties are 
then given by 

A{d,\) = ±l, (2.1) 
with measurement direction a on Alice's side, and 

B{b,\) = ±l, (2.2) 

— * 

with measurement direction b on Bob's side. Locality requires that the outcome of 
A does not depend on b and B does not depend on a. Without loss of generality we 
can say that the value A is achieved with probability density p(A) > obeying 

j piX)dX = 1 . (2.3) 

There is no quantum world. There is only an abstract physical description. It is wrong to think that the task 
of physics is to find out how nature is. Physics concerns what we can say about nature." - Niels Bohr 
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We define a correlation function for the joint spin measurement, that yields the value 
1 when the spins are parallel and —1 when the they are antiparallel 

(2.4) 



C yA{d, A), 5(6, A) j = A(a, A) • B{h, A) . 

The expectation value of this quantity therefore is 

E{d,b) = j p{\)A{d,\)Bib,X)dX . (2.5) 

If we combine different expectation values of different measurement directions in the 
following way 



E{d, h) - E{d, h') 



p(A) A{d,\)B{h,\) - A{d,\)B{h',\) dX 
p{X)A{d, X)B{b, A) [l ± A{d', X)B{b', A) 



- J p(A)A(a, A)5(6', A) 1 ± A{a', X)B{b, A) 
then the absolute value yields 



dX 



dX 



(2.6) 



\E(d,b) - E{d,b')\ < 



+ 



p(A) l±A{a',X)B{b',X) 
p{X) \l ± A{a', X)B{b, A) 



dX 



dX 



(2.7) 



= 2±\E{a',b') + E{a',b)\ . 

We obtain the CHSH inequality by rewriting this in the form 

\E{d, b) - E{d, 9) + E{a', 9) + E{9', b)\ < 2 . (2.8) 

The derivation reveals that every local realistic theory has to fulfill the CHSH in- 
equality. Quantum mechanics however predicts a violation under particular circum- 
stances. Consider the quantum mechanical expectation value of E{d, b) for the Bell 
state |vl;-) = -i=(|t|)-||t)) 

E{d,b) = {^-\d- a^b- al^^-) , (2.9) 

with unit vectors a and b. A short evaluation yields 

E{d,b) = -d-b= -cos{a- (3) , (2.10) 

wherein the angles a and /? substitute the vectors a and b. Hence, we can write the 
CHSH inequality in the form 

I - cos{a - (3) + cos{a - (3') - cos{a' - (3) - cos{a' - (3')\ , (2.11) 

wherein a, a', (3 and (3' are the angles of the four vectors d,a',b and b' in a plane. 
For angles obeying \a — l3\ = \a' — /5| = \a' — P'\ = j and \a — P'l = ^ we find 

I — cos{a — P) + cos{a — (3') — cos{a' — (3) — cos{a' — [3')\ 

= \-- — =2V2>2. 2.12 

'2222' ^ ' 

As a result we have found that any local realistic theory cannot reproduce the 

statistics of quantum mechanics. Experiments with entangled photons confirm the 

violation of the CHSH inequality (see [2], [22] )• 
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2.2.2 Horodecki violation criterion 



The example with the singlet state I^E'") demonstrates the violation of the CHSH 
inequality. Now we want to work out the whole set of states causing violations. For 
an arbitrary mixed state p the expectation value E{a, b) is 



E{a, b) = Tr{p a - a®b-a) 
Thus, we can write the Bell inequality in the form 



Tr ( p 



< 2 



a- a {b -U) ■ a + a' ■ a ® {b + b') ■ a 
We call the operator in square brackets the Bell operator B 

B{a,a',b,b') = a-a® {b - 9) ■ a + a' ■ a ® (b + 9) ■ a 



(2.13) 



(2.14) 



(2.15) 



It is obvious that for a given state p one has to find the right vectors a, a', b and 
b' to show possible violation. Only if the global maximum with respect to all Bell 
operators B{a,a' ,b,b') is less or equal to 2, then the inequality is preserved 



max 

d,a' ,b,b' 



Tr I p 



d-a®(b~9)-a + a'-a® (b + 9) ■ a 



< 2 . 



(2.16) 



Since the Bell operator B{d, a', b, b') is expressed in terms of Pauh matrices, it makes 
sense to express p in the way we introduced it in (11.161) 



1 / ' 

p=-|l(g)l+f-(T®l + l(g)S-a+ ^ tnm(^n ® (^n 
\ n,m=l 

The computation of \Tr{pB)\ then yields 

a ■ T{p){b- 9) + a' ■ T{p){b + 9) 

wherein T(p) denotes the 3x3 correlation matrix with the coefficients tr, 
replace b — V and 6 + 6' by mutually orthogonal unit vectors c and d 



(2.17) 

(2.18) 
We 



b-b' = 2smec' b + b' = 2cos^c . (2.19) 
Now the maximum of \Tr{pB)\ has to be determined with respect to 6, a, a', c and 



max 2 



sin6'a ■ T{p)c' + cos6a' ■ T{p)c 



(2.20) 



The scalar products are maximal for parallel vectors. Consequently, we choose the 

_ Tip)c' . -} _ T{p)c 



unit vectors Emax = and a'max 



sin^||T(p)c^|| + cose||T(p)^| 
Maximisation with respect to 6 yields 



max 2 



max2VI|T(p)c'P+ ||T(p)c1|2 . 

c,c' 



(2.21) 



(2.22) 
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Since c and d are mutually orthogonal unit vectors, the expressions ||T(p)c'|p = d • 
T'^{p)T[p)d and ||Tpc||^ = c- {p)T [p)c are maximal when d and c are eigenvectors 
of the two largest eigenvalues Ai(p) and A2(p) of K{p) = T'^{p)T{p). Thus the CHSH 
inequality reads 

2VAi(p) + A2(/)) <2 . (2.23) 

Now we have proven that if and only if the two largest eigenvalues of the matrix 
^ip) = T^ip)Tip) comply 

Ai(p) + A2(p) > 1 , (2.24) 

then the state p violates the CHSH inequality for particular vectors d,a',b and b'. 
This necessary and sufficient condition was found by the Horodeckis in 1995 [33] . 

2.3 A Bell inequality for bipartite qudit systems: CGLMP 

The CHSH inequality is a Bell inequality for bipartite qubit systems. In this section 
we present Bell inequalities for bipartite qudit systems based on logical constraints 
established by Collins, Gisin, Linden, Massar and Popescu (CGLMP) [31] . 

2.3.1 Derivation of the CGLMP inequality 

Once again, we consider the standard situation with the two parties Alice and Bob 
spatially separated, only having access onto their local Hilbert spaces = = 
of dimension d, while the composite Hilbert space is once again TC^^ = ® C^. 
Consequently, measurements on each side have d possible outcomes. Let us assume 
again that we could reproduce the statistical predictions with a local realistic theory. 
In analogy to the CHSH Bell experiment, each party has two apparatuses with 
different settings. Due to the fact that in a local realistic theory all observables 
have definite values simultaneously, the state of a system induces a probability 
distribution of form 

P{A,=j,A2 = k,Bi=l,B2 = m) (2.25) 

that apparatus Ai gives measurement result j G {0,..,d— 1}, apparatus A2 gives 
G {0, .., (i — 1} and so on. In sum, we have d'^ values determining the statistics of 
the system. As usual we normalise these probabilities 

J2 Pi^i = J, A2 = k,B,=l,B2 = m) = l . (2.26) 

jklm 

If we are only interested in measurement results of certain apparatuses then we 
have to sum over all ignored observables. For example the probability of Ai giving 
j and Si giving / is P{Ai = j, B^ = I) = J^km^i^i = ^2 = k,Bi = I, B^ = m). 
We introduce some variables r',s',t' and u' defined by relations of measurement 
outcomes 



r' 


= B, 


-A, 


= l~J , 


(2.27) 


s 


= A2 




= k-l , 


(2.28) 


t' 


= B2 


-A2 


= m — k , 


(2.29) 


u 


= A, 


~B2 


= j — m . 


(2.30) 
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The definition induces the constraint 

{r +s' + t' + u) modd = . (2.31) 

We introduce the quantity / defined by 

/ = P{Ai = Bi) + P(5i = ^2 + 1) + P{A2 = B2) + P{B2 = Ai) , (2.32) 

with the probabihties P{Aa = Bh + k) that the outcome of Aa differs from outcome 
Bb by k 

d-l 

P{A, = Bb + k) = J2 Pi^a = (j + k) mod d, Bb = j) . (2.33) 

j=0 

Then I is the sum of the probabilities that r' = 0,s' = —l,t' = and u' = 0. 
Nevertheless, because of the logical constraint r' + s' + t' + u' = only three of those 
four relations can be vahd. Thus we have 

J < 3 (2.34) 

for all local realistic theories. We introduce another quantity 

/s = + [P(Ai = 5i) + PiB^ = A2 + 1) + P(A2 = P2) + PiB2 = A^)] (2.35) 
- [P{A, = B,-1) + P{Bi = A2) + P{A2 = P2 - 1) + P{B2 = Ai- 1)] 

Consider the following table which shows how the logical constraint affects the max- 
imum of J3. 



+1 






1 








r' 


s' 


t' 


u' 


r' 


s' 


t' 


u' 









-1 








1 





1 


1 








♦ 


♦ 


♦ 


♦ 








^ h 


= 2 


♦ 




♦ 


♦ 




♦ 






^ /3 


= 2 


♦ 


♦ 




♦ 






♦ 




^ /3 


= 2 


♦ 


♦ 


♦ 










♦ 


^ /3 


= 2 



In consequence 

/3<2. 

Additional quantities Id can be introducecd 

f'^' 2k 

h= {^-jz^){ + [P{M = B^ + k) + P{B^ = A2 + k + l) 

+ P{A2 = B2 + k) + P{B2 = Ai + k)] 
- [P{A^ = B^-k-l) + P{Br = A2 - k) 
+ P{A2 = B2-k-l) + P{B2 = Ai-k 

For local realistic theories the upper bound is 2 

^the bracket [ ] stands for the floor function 



(2.36) 

(2.37) 

-1)]} 

(2.38) 
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r 


A 

= Ai 


- Bi , 


s 


= B, 


-A2-I 


t 


= A2 


-B2 , 


u 


= B2 





Proof: We introduce the variables 

(2.39) 
(2.40) 
(2.41) 
(2.42) 

that obey 

(r + s + t + u + l) mod d = . (2.43) 

Since all computations are done modulo d, without loss of generality we can restrict 
r,s,t and u to lie in the interval 

-[^]<r,s,t,u<[^-^]. (2.44) 

Id can be written as a function depending on r,s,t and u in the following way 

Id = f{r) + f{s) + fit) + f{u) , (2.45) 

where f{x) is given by 

-2x 



+ 1 a; > 

X <0 . 



/(^) = < Ux d+1 (2-46) 



d-1 d-1 

We now have to check all combinations of algebraic signs of r, s, t and u. With help 
of fl^:iH|) and (ICTj) we find 

• r, s,t,u > : 

^r + s + t + u + l = d^Id = 2 

• three of r, s, t, u are > and one is < : 
either r + s + t + u + l = d^Id = 2 
orr + s + t + M + l = 0^/d = j^ 

• two of r, s, t, u are > and two are < : 

• one of r, s, t, m is > and three are < : 
either r + s + t + u + I = Q ^ h = 
01 r + s + t + u + l = -d^ Id= ^ 



r, s,t,u < : 



r+s+t+u+l 



-2(d+l) 



^ -a d-1 

Hence, < 2. □ 
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In the following, we investigate the predictions of quantum mechanics for this family 
of inequalities. First, we have to insert the quantum mechanical probabilities 



(2.47) 



wherein and {lOsfc} denote orthonormal eigenvectors of the observables 

with fc, / G {0, .., — 1}. To confirm the violation of these inequalities consider the 
maximally entangled state 



^ d-l 



(2.48) 



s=0 



and orthonormal eigenvectors 



I Aa 



10 



Bb 



1 (2i,i , , 



A ' 



-1 + I3b) \s) 



'B ' 



(2.49) 
(2.50) 



with ai = , ^2 = I , /3i = ;i and P2 = Then 02.471) becomes 

=Tr{\no,o) {^o,o\ \k)^Jk\^^^\l)^, 

= mo,o\\k)Aa®\i)Bf 

'd-1 



1 
1 

1 



(^0,0 1 



d-1 



d-l 

E 

,r=0 



exp 



2m 



r{k + aa) \r) 



E 

,m=0 



exp 



27ri 



m(-/ + (3b) \m)^ 



s,r,m=Q 
d-l 



2m 



^ (Ma \^(Mb l"^exp ( —r{k + aa) ) exp ( —m{-l + /Sb 



2m 



s=0 



2m 



exp I -^s{k -l + aa + (3b 



_ sin^ {Tr{k - I + aa + (3b)) 
d? sin^ (7r(A; — / + + (3b/ d)) 

With the values of ai, 0:2, A and /32 given above we find 



P{Aa = k,Bb = l) 

Inserting this into / and Id leads to 
2 



1 



2d^ sin^ (7r(fc + + /^b) /d) 



rf2 



E 

k=0 



2k 



d-l) \ sin2(^(A; + i)) sin2(^(A; + |)) 



(2.51) 

(2.52) 

(2.53) 
(2.54) 
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We obtain the following values for / and Id 



d 


2 


3 


4 


5 


6 


7 




oo 


I 


3.41421 


3.31738 


3.28427 


3.26908 


3.26086 


3.25592 




3.24228 


violation [%] 


13.8071 


10.5793 


9.47559 


8.96922 


8.69533 


8.53059 




8.07593 


Id 


2.82843 


2.87293 


2.89624 


2.91054 


2.92020 


2.92716 




2.96981 


violation [%] 


41.4214 


43.6467 


44.8122 


45.5272 


46.0102 


46.358 




48.4906 



As we can see, the CGLMP inequalities / < 3 and Id ^ are violated. The 
table also reveals that by using Id instead of / we achieve much stronger viola- 
tions. This is in close relation to the resistance against noise. Before we analyse 
this, we present a concise notation for / and Id- Due to fl2.47p the probabilities 
P{Aa = Bb + k) become 



d-l 

Bb) 



P{A, = Bb + k) = Y,Tr (p |(j + k) mod d)^^ {{j + k) mod d\^^ ® {j\ 

j=0 



We rewrite / by exploiting linearity of the trace in the following way 
/ = P(Ai = 5i) + P{A2 = B2) + P(fii = A2 + 1) + P{B2 = Ai) 

(d-l 
P^{\j)Ai O'Ui ® \3)bi + \3)a2 {3\a2 ® \3)b2 01b2 
i=o 

+ |(j - 1) mod d)^2 ((j - 1) mod d\j^^ ® \i) (jI^^ + Ij)^^ (jI^^ ® \]) ^2 (jlij2 

= Tr{pBi) . (2.55) 

We identify the sum of operators as our new Bell operator Bj for /. In the same 
way we can define Bell operators Bj^ for all quantities Id so that 

Id = Tr{pBiJ . (2.56) 

Now, let us consider the state |f2o,o) in presence of uncolored noise 

p=(l-r)|fio,o)(^^o,o|+r^, (2.57) 

wherein r G [0, 1] is the amount of noise. We obtain 

/ =Tr (pBi) = (1 - r)Tr (|^]o,o) (^^o,o| ^i) + ^TrBj , (2.58) 

Id =Tr ipBj,) = (1 - r)Tr (|fio,o) (^^o,o| + ^TrBj, . (2.59) 

Bi contains 4d projectors with prefactor 1, hence TrBi = 4d. Whereas the addends 
in Bi^ contain Ad projectors with prefactor (l — and Ad projectors with prefactor 
— (1 — hence TrBj^ = 0. In consequence, we obtain 

I =(1 - r)Tr {\Qo.o) {^o,o\ Bi) + r- , (2-60) 

d 

Id={^-r)Tr{\no,o){no,o\Bj,) . (2.61) 
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Above a certain amount of noise r = Vmax, the state p ceases to violate the Bell in- 
equality. The values are given by I{rmax) = 3 and Id{,fmax) = 2, and are summarised 
in the following table. 



d 


2 


3 


4 


5 


6 


7 




oo 


I • ^max 


%] 


29.2893 


15.9965 


12.4446 


10.8979 


10.0555 


9.5332 




7.47246 


Id • ^max 


[%] 


29.2893 


30.3848 


30.9450 


31.2843 


31.5116 


31.6744 




32.6557 



We conclude that one should give preference to the inequalities Id ^ due to 
their higher resistance to noise. 



2.3.2 Optimisation of the Bell operator 

The orthonormal bases given in fl2.49p are optimal for the state |fio,o)- In other 
words, there exists no better Bell operator for |f2o,o) causing higher violation. This 
has not been proven so far, though numerical optimisation indicates that this is true 
(see [3l],[35]). The lack of a proof lies in the fact that the finding of an optimal 
Bell operator for an arbitrary p is a nonlinear optimisation problem: The finding of 
solutions for grad(Tr(pi3)) = with respect to all vector components is a nonlinear 
task and is further complicated by nonlinear constraints since the vectors {1/^)^^} 
and {lOfib} have to form an orthonormal basis (observables are hermitian operators 
and their eigenvectors are orthogonal) 

{k\Aa \j)Aa = 4, , (2.62) 
mBb\^)Bb = Sim ■ (2.63) 

Using the Weyl operators or the Gell-Mann matrices to express p and B does not 
lead to a simplification like in §2.2.21 for the CHSH inequality. For this reason, we 
have developed a numerical optimisation algorithm which reliably finds optimal Bell 
operators for any p. Let us first investigate the number of variables of the Bell opera- 
tor B. We have 4 orthonormal bases , {1^)^21 ' {IObiI {11)^2} with each 
d vectors. Each vector is an element of and can be described by 2d real numbers 
(d vector components, each described by 2 real numbers, one for the amplitude and 
one for the phase). Altogether we have 4 x d x 2d = Sd"^ real variables. Orthonor- 
mality restricts the values of the variables, thus we have to optimise B under the 
constraints (12.621) and (12.631) . For all practical purposes however this way is quite 
impractical. We now show how to satisfy the constraints by choosing the right set 
of variables. We begin with some mathematical considerations. Any basis trans- 
formation can be realised via a unitary transformation, i.e. any basis {\k)} of 
can be transformed into any other basis of by means of a certain unitary 
transformation {]/;;')} = {U\k)}. Since it is our goal to find the right orthonor- 
mal bases , {1^)^12} ' {IObiI {I0_b2}' "^^^ select arbitrary orthonormal 
bases and seek the unitary transformations Uai,Ua2,Ubi and Ubi that maximise 
Id and /. These transformations are elements of the unitary group U{d). Due to 
the fact that global phases do not affect probabilities and that U{d) can be written 
as a semidirect product U{d) = SU{d) x f/(l), it suffices to regard the special uni- 
tary group SU{d). We use the generalised Euler angle parametrisation of this 
group, wherein the Euler angles embody all degrees of freedom (see [36j,[37j). For 
this kind of parametrisation the following antisymmetric and diagonal generalised 
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Gell-Mann matrices (11. 38111. 39p are needed 



'k-2 



A 



fc2-l 



A 



fc-2 



{k-l)k 



Y^\j){j\-{k-i)\k-i){k-i\ 



0=0 



A 



fc2 + l 



= A' 



Ok 



-z|0) {k\^i\k) 



2<k<d 



l<k<d-l . 



In this notation the exphcit form of the parametrisation oiU E SU{d) by generahsed 
Euler angles {ttj} reads (see v^t^ 



U 



'd-2 /d-x 



JJ exp(iA„2_ia^2_(rf+i_„)) 



n=2 



lx=0 \k=2 

with A{k,j{x)) and j(a;) given by 

MkJix)) = exp(iA3a(2fc-3)+i(x)) ■ exp(iA(fc_i)2+ia2(fe-i)+i(m)) , 



(2.64) 



(2.65) 







x = 
X > 



1=0 



For instance, for = 2, 3, 4 we find 



U e SU{2) 
U e 5f/(3) 

U e SU{A) 



U =exp(zA3ai) 
U =exp(zA3ai) 
•exp(zA3a5) 
U = exp(zA3ai) 
•exp(zA3a5) 
■exp(zA3a9) 



exp(iA2a2) ■ exp(2A3a3) , 
exp(zA2«2) ■ exp(zA3a3) ■ exp(zA5a4) 
exp(zA2a6) ■ exp(zA3a7) ■ exp(zA8a8) , 
exp(iA2a2) ■ exp(2A3a3) ■ exp(iA5a4) 
exp(zAio«6) ■ exp(zA3a;7) ■ exp(zA2a8) 
exp(zA5aio) ■ exp(iA3Q;ii) ■ exp(iA2a;i2; 



(2.66) 



(2.67) 
(2.68) 



• exp(2A3ai3) ■ exp(iA8ai4) ■ exp(iAi5ai5) . (2.69) 

The parametrisation guarantees comphance with the constraints (12.621) and (12.631) 
and reduces the number of variables to 4((i^ — 1) without discarding any solution. 
However, the optimisation of B with respect to the Euler angles still requires a nu- 
merical optimisation algorithm. While popular algorithms like Differential Evo- 
lution, Simulated Annealing or gradient based methods have been very time- 
consuming, the Nelder-Mead method [38j has performed this task relatively fast 
and reliable. In general this method can be advantageous when the number of 
variables is very large, because it manages to find a maximum without computing 
derivatives, which can possibly be computationally intensive. For maximisation the 
algorithm proceeds as follows: Assume a function f{xi, .., x„) which has to be max- 
imised with respect to n variables y = (xi, ..,x„). At the beginning a simplex with 
n + 1 vertices yi, ..,yn+i is created at random. In the first step (1) the vertices are 
being arranged by their values f{y) and labeled according to 



f{yi) > f{y2) >■■■> /(z/n+i) . 



(2.70) 



°Note that our notation differs slightly. We chose an indexing where all indices are increasing and the sequence 
of the product is fli^i — ' ^2 ■ ■ ■ ^jv 
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According to this, f{yi) is the best and /(j/„+i) the worst vertex. Then a reflection 
of the worst point through the centroid of the remaining n points ?/o = ~ Yl^=i Hi is 
performed 

y^=yQ + a{yo - yn+i) , (2.71) 

where a > is called the reflection parameter. The next step of the procedure 
depends on the value f{yr)- 

• If f{yr) is better than but not better than f{yi), i.e. f{yi) > f{yr) > 
f{yn) then a new simplex with |/„+i substituted by yr is built and step (1) is 
repeated. 

• If f{yr) is the best point, i.e. f{yr) > f{yi) then an expansion ye = yo + jiyo — 
yn+i) with expansion parameter 7 > is performed. If f{ye) > fijir) then 
a new simplex with yn+i substituted by ye is built and step (1) is repeated. 
Otherwise yr is used instead of ye- 

• If does not result in an improvement, i.e. fiyr) < /(yn+i) then a contrac- 
tion ye = yn+1 + p{yo — yn+i) with contraction parameter p > is performed. 
If f{yc) > /(l/n+i), then a new simplex with yn+i substituted by yc is built and 
step (1) is repeated. Else the simplex is being shrinked: The vertices are being 
replaced by y[ = yi and y'i = yi + cr{yi — yi) with shrink parameter cr > for 
all i G {2, ..,n + 1}, afterwards step (1) is repeated. 

The rules are repeated until the convergence criteria \f{y'i) — f{yi)\ < Ci and 
Wy'i ~ yi\\ < C2 are satisfied, where y[ is the new and yi the old best point, and 
Ci,C2 > are constants depending on the desired precision. As we can infer from 
the rules, the Nelder-Mead method is a hill climbing algorithm, which has the dis- 
advantage that it could converge at a local maximum. We cannot completely avoid 
this problem; however, by varying the starting points we will find the global maxi- 
mum in all likelihood. 

The optimisation procedures for the Bell operators iS/g and Bj^ via the Euler angle 
parametrisation and the Nelder-Mead method have been realised in MATHEMAT- 
ICA 6 and can be found in the appendix [B] and O To achieve good results it was 
necessary to find proper values for the parameters a,7,p and a. With the values 
a = 1.6,7 = 1-6, p = 0.8 and a = 0.8 an agreeable compromise between robustness 
against local maxima and time exposure was attained (with the standard values the 
algorithm converges at local maxima more often). We have chosen to execute the 
algorithm ten times with different starting simplices to guarantee that the global 
maximum is obtained. The accuracy/precision goal has been set to MATHEMAT- 
ICA 6 standard values Ci = C2 = 10~^. 

2.3.3 Properties 

Before we discuss some properties of the CGLMP inequalities we show that the 
inequality < 2 is equivalent to the CHSH inequality when the regarded Hilbert 
space is two dimensional TC^ = Ti.^ = C^. Correspondingly, all properties of the 
CGLMP inequalities Id ^ likewise hold for the CHSH inequality. 
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The proof of equivalence is straight forward. We write all terms of I2 explicitly 



+ [P{Ai = fii) + P{B^ = A2 + 1) + P(A2 = B2) + P(fi2 = ^1)] 

- [P{Ai =Bi-l) + P{Bi = A2) + P{A2 = B2-1) + P{B2 = 
+ [P{A, = 0, 5i = 0) + P{Ai = 1, 5i = 1) + P{A2 = 0, 5i = 1) 

+ P{A2 = l,Bi = 0) + P{A2 = 0, 52 = 0) + P{A2 = 1, ^2 = 1) 
+ P{Ai = 0, 52 = 0) + P{Ai = 1,B2 = 1)] 

- [P{A, = 0, 5i = 1) + P{Ai = 1, 5i = 0) + P{A2 = 0,Bi = 0) 
+ P{A2 = 1, 5i = 1) + P{A2 = 0, 52 = 1) + P{A2 = 1, 52 = 0) 
+ P(Ai = 0, 52 = 1) + P{A, = 1, 52 = 0)] . 



1)] 



Since the expectation value fl2.5l) can be written as E(A, B) = f p{X)A{\)B{X)dX = 
P{A = 0, 5 = 0) + P{A =1,5 = 1)- P{A = 0, 5 = 1) - P{A = 1, 5 = 0) we can 
summarise the terms of I2 and find 



I2 = E{A,, 5i) - 5(^2, 5i) + 5(^2, 52) + E{A^, B2) < 2 . (2.72) 

Due to the symmetry of positive and negative terms of I2 for dimension two, the 
inequality I/2I < 2 holds too. Hence, this is exactly the CHSH inequality with A 
and 5 interchanged (which makes no difference). 

Now we quote some important properties: The CGLMP inequalities have been 
proven to be tight Bell inequalities (see |39]). For the understanding of this, we 
have to return to the probability distribution fl2.25p on which the CGLMP inequal- 
ities are based on. As noted before, there are d'^ values determining the statistics 
of a local realistic theory. In contrast, in quantum mechanics we have probability 
distributions P{Aa = k, Bf, = I), with k,l E {0, ..,d — 1} for each of the four settings 
{AiBi), {A1B2), {A2B1) and {A2B2). As we have seen, those cannot be reproduced 
by summing over all unregarded observables of P{Ai = j, A2 = k,Bi = I, B2 = m). 
To conclude, we only have d"^ probabilities for each setting and Ad"^ in total. The 
quantum mechanical probability distributions are restricted by normalisation 

d-l 

Y^P{Aa = k,Bb = l) = l a, 6 =1,2, (2.73) 

k,l=0 

and the non-signaling condition (quantum mechanics cannot be used for superlumi- 
nal communication) 

d-l d-l 

J2p{Aa = k,B, = l) = Y,PiAa = k,B2 = l) a = 1,2 (2.74) 
1=0 1=0 



(it is clear, that the same has to hold for A and B interchanged) . 
What have we gained? The quantum mechanical probabilities can be specified by Ad"^ 
values, which can be regarded as entries of a vector P G 7?.'^°'^ . Due to the constraints, 
all physical relevant vectors lie in an affine space of dimension Ad'^ — Ad (see P9]). 
In this context, the normalisation coming from the local realistic description of the 
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system restricts a vector P to belong to the convex hull of d'^ vectors Gj G TZ^'^^ 

d-l 

P{Ai=j,A2 = k,Bi = l,B2 = m) = l, (2.75) 

j,k,k,m=0 

^Plr = Y1 c^G. q > and ^ q = 1 . (2.76) 

The extremal points Gj G TZ'^'^^ then are the generators of a convex polytope. Such 
a polytope can also be described by facets and their half-spaces. Facets induce 
inequalities of form 

P G conv {{Gi}) ^ X, • P < Vi G {1, 2, .., n} , (2.77) 

where n is the necessary number of facets to define the polytope. Hence, with this 
set of inequalities we can determine all states that are in contradiction to local re- 
alism in an experiment with two apparatuses on each side. The problem involved 
is, that computing the facets of a high-dimensional polytope is a very difficult task 
that has only been completely solved for d=2 (see [10] )• Let us now compare this 
approach with our well-known Bell inequalities. We know that all states that vi- 
olate a Bell inequality do not belong to the convex hull of {Gj}. However, states 
that fulfill a Bell inequality do not necessarily belong to it. This means that Bell 
inequalities define a region (half-space, sphere, polytope, etc.) that contains the 
convex polytope. The term tight Bell inequality is used if the boundary of the 
generated region at least partially coincides with at least one of the facets of the 
polytope. Hence, the inequalities 02.771) are themselves tight Bell inequalities as well 
as the CGLMP inequality given by Jg^. For the CGLMP it has been shown in [30] 
that it coincides with a family of equivalent facets. However, we do not know if 
the CGLMP inequality coincides with all facets of the polytope. Even though this 
might be the case, we cannot infer from this that all non-local states can be found 
with help of the CGLMP inequality because increasing the number of observables 
per party could define an improved polytope that could enable us to find even more 
non-local states. Please refer to [41j for a more detailed discussion. 
Now let us investigate another property of the CGLMP inequality concerning the 
maximal violation. As we stated before, for the maximally entangled state |fio,o) 
the violation is maximal when the measurement setup is configured according to 
(12.491) . However, in this case the largest eigenvalue of the corresponding Bell opera- 
tor Bi^ is larger than the value Id (|f2o,o)) for all d > 2. The largest eigenvalues are 
summarised in this table (see 1^): 



d 


2 


3 


4 


5 


6 


7 




2.8284 


2.9149 


2.9727 


3.0157 


3.0497 


3.0776 




2.8284 


2.8729 


2.8962 


2.9105 


2.9202 


2.9272 


difference [%] 





1.4591 


2.6398 


3.6133 


4.4345 


5.1411 



wherein {"^mv) denotes the eigenvector of Bj^ with the largest eigenvalue. Since all 
maximally entangled states are equivalent in terms of local unitaries (as discussed in 
§1.1. 3p the state |\E'm^) must be non-maximally entang leci. For instance, for qutrits 

^mii stands for "maximal violation" 

IHo^o) = Ua Cgi Ub I'^mv) then we could use the Bell operator Ua Cs) UbBi^U\ (gi Ug to obtain (|f2o,o)) = 
Id il^mv)) in contradiction to obtained values. 
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the eigenvector is the non-maximally entangled state I'^mv) = (|00) + 7 |11) + |22)) 

with n = 2 + 7^ and 7 = (vTT— v^) /2. Knowing that the states I'ifmv) yield 
stronger violations, we could ask if the above given values are the highest obtain- 
able for Jrf because the measurement configuration fl2.49p has only been intended to 
be the best for |fio,o)- Numerical investigations with varying states and measure- 
ment configurations show that these values are indeed the highest (see |13]). How 
can we interpret this result? Does the higher violation and the consequential higher 
resistance against noise signify that the states {"^mv) are more non-local than |f2o,o)? 
This and further questions on the comparison of nonlocality and entanglement will 
be discussed in the next section. 

2.4 Remarks on nonlocality and entanglement 

A state p is said to be non-local if it violates a Bell inequality. Such a state 
is necessarily entangled, which follows from the fact that separable states cannot 
violate any Bell inequality. This can easily be seen by generalising the probability 
distribution (12.251) . Recall that ideal Bell inequalities correspond to the facets of the 
polytope of local realistic correlations (12.771) . The number of observables on each 
site depends on the regarded Bell inequality but we do not want to restrict ourself 
to the case of two observables. Thus, we generalise (I2.25P to 

P(Ai = J, ..,Ar = k, Bi = l,..,B, = m) , (2.78) 

with r,s E N. According to (12.761) . the normalisation defines a polytope and all 
local states belong to it. It is not difficult to recognise that any separable pure state 
belongs to the polytope: Since all probabilities are of product form P{An = a, B„i = 
b) = P{An = a) ■ P{Bm = h) (which is not the case when a state is entangled) we 
can easily construct the above stated probability distribution 

P(Ai = J, .., A, = /c, 5i = /,.., 5, = m) 
=P{Ai =])■■■ P{Ar = k)P{Bi = /)■■• P{Bs = m) , (2.79) 

which of course obeys the normalisation relation due to normalised state vectors. 
Since separable mixed states are only convex combinations of pure ones we con- 
clude that they belong to the polytope as well. This completes the proof that 
separable states cannot violate any Bell inequality. As a consequence, the operator 
[2 ■ 1 — Bj^] G is a entanglement witness because Tr{pBiJ < 2 holds for all 

separable states. However, in general this is not an optimal witness, even if Bi^ is 
optimised for p. 

Let us go back to the main subject of this section. Entanglement is a necessary 
component for nonlocality, however not all entangled states violate a Bell inequal- 
ity. This will become more obvious when we discuss the geometry of entangle- 
ment and nonlocality in ^ At this point, we mention one famous example of an 
entangled state which does not violate the CHSH inequality: The Werner state 
pw = (^1 + (for d=2 equivalent to the isotropic state) is entangled 

for ^ < p < 1, however violation of the CHSH inequality is obtained only for 
< p < lH. This means that for | < p < ^ the statistical predictions can be 
reproduced by a local realistic theory, even though the state is entangled. This 
brings us to the notion of hidden nonlocality. We can claim that any distillable 

^The values of the parameter p can easily be verified with the PPT criterion and the Horodecki violation criterion 
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state ( §1.4.3p contains some amount of nonlocality since we could perform local op- 
erations and classical communication ( §1.4.2p until nonlocality is revealed. Despite 
this argument, the distinction between nonlocality and entanglement is not entirely 
clarified. Besides the fact that we do not know if any entangled state behaves non- 
locally in a particular measurement situation, there are some well-known examples 
that indicate that they might not be exactly the same. Eberhard has shown that 
non-maximally entangled states require lower detection efficiencies than maximally 
entangled ones, in order to close the detection loophole (see [12] )• We have also 
seen that some non-maximally entangled states seem to be more non-local than the 
Bell states, due to their higher violation of the CGLMP inequality. With regard to 
this issue, one might argue that the resistance against noise is not a good measure 
of nonlocality. Some remarks on this can be found in a publication by Acin, Durt, 
Gisin and Latorre [13]. Nevertheless, they also suggest that a Bell inequality with 
more than two measurement instruments on each site could avoid this peculiarity. 
Another hint that entanglement and nonlocality are different resources has been 
found by Brunner, Gisin and Scarani (see 01]). They have shown that the simu- 
lation of non-maximally entangled states via non-local machines requires more 
resources than the simulation of a Bell state. The hypothetical non-local machinJ^ 
was constructed to obtain the algebraic bound of the CHSH inequality without vio- 
lating the non-signaling condition. Alice and Bob each have an input (x and y) and 
an output (a and h). Alice can choose the value of a; G {0, 1} and gets a G {0, 1} in 
return, while Bob can choose the value of y G {0, 1} and gets 6 G {0, 1} in return. 



Alice 



X e {0, 1} 



a e {0,1} 



non-local machine 

(a + b) mod 2 — x ■ y 



ye {0,1} 



be {0,1} 



Bob 



Figure 2.2: Schematic illustration of the non-local machine 

The machine creates a random output P{a = 0) = P{a = 1) = P{b = 0) = P{b = 
1) = though with a and b correlated according to the rule 

(a + b) mod 2 = x-y . (2.80) 

For instance, if x = y = then the output of the machine is a = 6 = in half the 
cases and a = 6 = 1 in the other half. Due to the local randomness this box cannot 
be used for signaling. Let us investigate the violation of the CHSH inequality. We 
assume that for any observable A{a) there is a corresponding value for x{a) G {0, 1}, 
and the same applies to B{b) and y{b) G {0, 1}. As before, we define a correlation 
function that yields the value 1 when the outputs a and b are identical, and — 1 
when they are different 

c{A(a),Bii))^[_\ (2.81) 

For the case A{a) x{a) = l,A{a') x{a') = 0,B{b) y{b) = and B{b') 
y{b') = 1 the expectation values of the correlation function are E{a, b) = E{a',b) = 

other common names are PR box (named after Popescu and Rohrlich) and non-local box 
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E(a',l/) = 1 and E{a,b') = —1. Under such circumstances, the CHSH inequahty 
yields the algebraic bound 

\E{a,b) - E{a,b') + E{a',b') + E{a',b)\ = 4 . (2.82) 

This is the most non-local behaviour of a system possible. Thus, from a informa- 
tion theoretical point of view, this machine can be seen as a resource of nonlocality. 
Inspired by this idea, Cerf, Gisin, Massar and Popescu have shown that quantum 
entanglement can be simulated without communication by use of the non-local ma- 
chine (see |15]). Subsequent work demonstrated that for the simulation of a Bell state 
only one non-local machine is required, while at least two are required to simulate 
non-maximally entangled states (see [Bj). This is another fact that strengthens the 
conjecture that there is a difference between nonlocality and entanglement. Further 
open questions concerning Bell inequalities, including also the experimental point 
of view, can be found in a publication by N. Gisin ^46j. 
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3 Geometrical aspects of bipartite systems 



As the state space in quantum physics is a complex Hilbert-Schmidt space 'M-Hs it is 
in general difficult to get a feeling for the properties of density matrices. However, in 
some cases they can be represented by vectors in a real vector space. For example, 
in §1.11 we have shown that density matrices can be described by real Bloch vectors. 
Due to the fact that the dimension of this vector is c?^ — 1, this representation is not 
very helpful if we intend to investigate bipartite systems. This is because even for 
the simplest case of two qubits the vector space is 15 dimensional. In this section we 
show that it is still possible to find attractive illustrations when it comes to studying 
entanglement and nonlocality of bipartite systems. 

3.1 Geometry of the two-qubit system - the tetrahedron 
3.1.1 Introduction 

As usual, we begin with the bipartite qubit system (the following approach can be 
found in [1^ and |1H])- According to fll.l6p . any density matrix acting on Ti.^^ = 
® can be written in the form 

1 ^ 

P=-(l(g)l+f-(7(8)l + l(g)S-a+ ^ tnm(^n ® (Tm) ■ (3.1) 

n,m=l 

Separability and nonlocality are invariant under local unitary transformations Ua ® 
Ub, therefore we can define equivalence classes of states 

[p] = {p\p =Ua<^ UbpU\ ® Ui) (3.2) 

that have the same properties concerning separability and nonlocality. We exploit 
the group isomorphism between SU{2) and SO (3) which induces that for any O G 
5*0(3) there exists a. U E SU{2) that obeys 

Un ■ aU^ = (On) ■ a . (3.3) 

Under a transformation Ua®UbpU\®U'^^ the vectors r, s and the matrix T = {tnm) 
then become 

r' =OAr , (3.4) 
s' =Oas , (3.5) 
r =0aT01 . (3.6) 

According to the singular value decomposition, there exist orthogonal matri- 
ces Oa and Ob so that T' becomes diagonal with real entries. We choose den- 
sity matrices p with diagonal T = diag{tii,t22,Hz) as the representatives of the 
equivalence classes [p]. Hence, for determining separability and nonlocality it suf- 
fices to investigate the set of states with diagonal T which we denote by V. We 
write the diagonal entries in a vector t = (^11,^22,^33), in this way a density ma- 
trix p G P is described by three vectors r,s,t G M^. Now, let us look through 
the elements of T>. Non-negativity is a necessary condition for a density ma- 
trix p, i.e. Tr d"^) (^|p) > for all |\E') G Ti.^^ . Consider the four projectors 
= |^+) (^+1 ,P2 = 1^") (*~| , A = 1'^'+) ('^'+1 and P4 = 1$-) ($-| given by the 
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Bell states = ^ (|01) ± |10)) and = ^ (|00) ± |11)). The Pauli matrix 
decomposition fll.23p yields rj = and Sj = for i G {1, 2, 3} and 

t1 = (+l,+l,-l) , (3.7) 

t2 = (-1,-1,-1) , (3.8) 

r3 = (+l,-l,+l) , (3.9) 

t4 = (-1,+1,+1) . (3.10) 

For two states p and p' given in the Pauli matrix decomposition we can compute 
Tr{pp') easily, because it simplifies to Tr{pp') = ^(1 + r ■ + s ■ s' + Tr{TT'^)). 
Thus, the four inequalities Tr (PnP) > with n G {1, ..,4} give 

1+^11+^22-^33 > , (3.11) 

1 - til - ^22 - ^33 > , (3.12) 

l+tll-t22 + t33>0 , (3.13) 

1-^11+^22 + ^33 >0 . (3.14) 

This restricts any t of p G P to belong to a regular tetrahedron spanned by 
the Bell states at the vertices (see Fig J3.1l) . Note that only in the case where 
r = and s* = the density matrices p G P are diagonal in the Bell basis 
{|^+) , 1^") , 1$+) , 1$")} and only then the constraints fl3?TTl) -f l3?TIl) are also suffi- 
cient for non- negativity. States with this property define a subset LMM C T) and 
are called locally maximally mixed because the respective reduced density ma- 
trices are maximally disordered Tr^^p) = TrB{p) = |l. The subset LMM C V 
contains only four pure states, namely the four Bell states, which can easily be seen 
by testing the purity condition 

Tr(p2) = i(l + 4+4 + 4) = l (3.15) 

^4+4+4 = 3 . (3.16) 

Geometrically spoken this defines a sphere with radius that intersects the tetra- 
hedron at the vertices, which therefore are the only pure sates of LMM. 



Pi = 1 








P2=\ 




' {'^ \ 




P3 = 


1$+) 


>($+ 




Pa = 


1$-; 


> 
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Figure 3.1: Illustration of the regular tetrahedron (green) 



3.1.2 Geometry of separable and entangled states 

Now, let us study the separability of p G LMM. In §1.2.21 we stated that the PPT 
criterion is a necessary and sufficient criterion for separability for density matrices 
acting on Ti^^ = C^®C^. In the Pauli matrix decomposition it is trivial to compute 
the partial transposition, ai is a symmetric matrix and therefore invariant under 
transposition, the same applies to which is diagonal. Only (T2 changes its algebraic 
sign (tJ = — (J2, since it is anti-symmetric. Consequently, partial transposition results 
in a reflection of the tetrahedron where the coordinates are transformed according 
to (til, ^22, ^33) (^11,-^22,^33)- All operators p G LMM belonging to the original 
tetrahedron are non-negative, while the reflected tetrahedron contains operators p 
with negative eigenvalues. This means that all states p belonging to the cross section 
of the two tetrahedron are separable, while the other states are entangled. This 
intersection is an octahedron (see Fig |3.2p . The extremal separable points have 
the coordinates ti/2 = (0, 0, ±1), ^3/4 = (0, ±1,0) and ts/g = (0,0, ±1) and are 1:1 
mixtures of two Bell states. A possible decomposition into a convex combination of 
separable density matrices can easily be found. Consider the mixture p = |(pi + P2) 

of the separable states pi = |(l + a-a)®|(l + 6-a) and p2 = |(1 — a-c?)®|(l — 6-a). 
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A short computation gives 

11- 
P = l^{pi + P2) = ^{t + a-a®h-a) . (3.17) 

With the vectors a = Cn and h = ±e„, where e„ stands for the unit vectors 
6*1 = (1,0,0), 62 = (0,1,0) and 63 = (0,0,1), we can exphcitly demonstrate the 
separabihty of the extremal points of the octahedron, thus all states within the 
octahedron because they are just convex combinations of these points. 

1$-) ($-| 




Figure 3.2: Illustration of the tetrahedron (green) and the octahedron of separable states (blue) 



3.1.3 Geometry of non-local states 

We intend to find the non-local states within the tetrahedron, i.e. states that violate 
the CHSH inequalit}l"1. For the bipartite qubit system this is not demanding because 
we can exploit the Horodecki violation criterion fl2.24p . First, we have to compute 
the eigenvalues of the matrix K = T^T. Since T is already diagonal we get Ai = 
^11; A2 = ^22 -^3 = ^33- ^ sufficient condition for nonlocality is that the sum of 
two eigenvalues is larger than 1 which implies the validity of at least one of the three 

"CHSH is equivalent to CGLMP for d=2 
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inequalities 



4+4>l, (3.18) 
4+4>i, (3.19) 
4 + 4 >i. (3.20) 

Each of the inequahties defines a cyhnder with radius 1 in the geometric picture and 
a density matrix p that hes outside of one of them is non-local. The union of the 
exterior regions of these three cylinders is the set of non-local states (see Fig J3.3l) . 

1$-) ($-| 




Figure 3.3: Illustration of the tetrahedron and the borders of nonlocality. States p beyond the 
red/meshed surfaces violate the CHSH inequality. 

This descriptively illustrates that not all entangled states violate the CHSH inequal- 
ity and that Bell operators Bchsh are non-optimal entanglement witnesses for LMM 
states. 
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3.2 Geometry of bipartite qudit systems - the magic simplex 
3.2.1 Introduction 

The derivation of the octahedron was direct and straight forward. In seeking a 
similar geometric object for qudits we definitely try to proceed analogously. The 
first step is to express p G 'Hjj^ in the generalised form of (11.161) which is fll.43p 



P 



\ i=l 



(3.21) 



1=1 



As we require a real vector space for a geometric structure, we use the generalised 
Gell-Mann matrices for {Fj}. In order to use the same notation as for qubits, we 
express p with generalised Bloch vectors ta, s\ and a correlation matrix T = {tij) 



l(g)l + fA-A(g)l + l 



SA ■ A + ^ UjVi ® 



(3.22) 



The next step would be the diagonalisation of the matrix T = (tij) via two orthog- 
onal matrices Oa and Ob- However, here we are confronted with the problem that 
the relation Uf^ ■ AW = {Of a) ■ A is invalid for Gell-Mann matrices. For example, 
consider the state p = |(1 + 6a ■ A) with = (0, 0, 1, 0, 0, 0, 0, 0) 



p= -(I + A3) 



f 

i 



(3.23) 



which is obviously a valid density matrix. Let us assume that there exists a unitary 
transformation p' = UpW = |(1 + Ub\ ■ AW) = |(1 + (OFa) ■ A) that corresponds 
to the orthogonal matrix O 








Vi/ 

This transformation turns the state p into p' = ■^{1 + b'j^ ■ A) 
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(3.24) 



(3.25) 



As we can see the eigenvalues of p and p' are not the same, which in other words 
means that the orthogonal transformation O results in a non-unitary transforma- 



tioio of p. Hence, the relation Uf^ ■ AW = (Of a) ■ A cannot be valid. Even if 



^Eigenvalues are invariant under unitary transformations 
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T can be diagonalised by product unitary transformations Ua ® Ub for a subset of 
LMM states there is another problem. Consider a focally maximally mixed state 
p e LMM El with diagonal T 

P = jj^i®i + XJ ^-r, ® r}j . (3.26) 

In the qubit case the eigenvectors of p are the Bell states for an arbitrary choice of 
{tii}, thus the non- negativity condition leads to a simple geometric figure, namely 
the tetrahedron. For qudits the situation is different because the eigenvectors of p 
depend on {ta} implying that the geometric figure given by non-negativity is much 
more difficult to determine. Due to this and the fact that the set of states whose 
correlation matrices T can be diagonalised via local unitaries are merely a subset of 
LMM, we choose another subset of LMM whose properties concerning entanglement 
are interesting and whose non-negativity is easier to handle. 

At this point we adopt the way of B. Baumgartner, B.C. Hiesmayr and H. Narnhofer, 
who introduced a generalisation of the tetrahedron for qudits |19]. As the octahedron 
is spanned by the four Bell states, the expansion for higher dimensional systems is 
the so called magic simplex spanned by generalised Bell states Pkj = \^k,i) {^k,i\ 
with \ilk,i) = 0^k,i ® 1) 1^0,0) given by the Weyl operators flL3ip 

d-l d-l 

W = { 5^ Ck,iPk,i I Ck,i > 0, 5^ Ck,i = 1} . (3.27) 

k,l=0 k,l=0 

This simplex is a convex set located in a c?^ — 1 dimensional hyperplane in a 
dimensional real vector space of hermitian operators spanned by the operators {Pk,i}- 
When the element p G W C TY^f is expressed in the basis ;)} G it becomes 
obvious that the only pure states in W lie at the vertices {Pk,i} (diagonal p). Before 
we investigate the states within the magic simplex let us specify the subset of LMM 
states belonging to it. One requirement for the existence of a representative in W 
surely is that p G LMM must be decomposable into orthogonal Bell states. However, 
we now show why this is only a necessary and not a sufficient condition. Consider 
the unitary transformation U = Ua® Ub = Xls k') (-^l ® 1 transforming a Bell state 
1^°'°^ ~ Vd ^« 1'^^ 1'^^ another orthonormal Bell state |$) = J2s \^') ® k)- 
Orthogonality implies TtUa = meaning the sum of all eigenvalues Ua has to be 
zero. Without loss of generality, we can set one of the eigenvalues 1 because we 
are free in choosing a global phase. This implies the rest of the eigenvalues for 
qubits and qutrits, while for d> A there are various ways for complying TtUa = 0. 
The Weyl operators W^^i have the eigenvalues e^'^'^^l'^ with 6 G {0, .., — 1} and if a 
state is decomposable into certain Bell states their intertwiners must have the same 
eigenvalues in order to have a representative in W. 



3.2.2 Symmetries and equivalences inside W 

We focus on local ( ant i-) unitary transformations U a ® Ub mapping W onto it- 
self. In particular we are interested in equivalence classes [p] = {p' G VV|p' = 
U A ® UbpU\ ®U'^b\ states within the magic simplex having the same properties 
concerning separability and nonlocality. These equivalences can best be studied 

13LMM = {p e 'HjPs\pA = TtbP = Pb= TtaP = \l} 
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when the group structure of the Weyl operators and the concept of a finite discrete 
classical phase space are used. The Weyl operators originate from the quantiza- 
tion of classical kinematics where they are used for translations between discrete 
states in a phase space. Each point within this phase space corresponds to an index 
pair {k, I), where Z G {0, .., o? — 1} denotes the quantum number of the position and 
k e {0, .., d — 1} the quantum number of the momentum. 




Figure 3.4: Illustration of the finite discrete classical phase space for rf = 3. 

Due to the fact that the Weyl operators have been used as the intertwiners between 
the vertices {Pk,i} of the simplex we can exploit this phase space structure for our 
investigations, i.e. to each point {k, I) of the lattice in the phase space we assign a 
Bell state We know by construction that any Bell state Pk^i can be mapped onto 
any Pfc/j/ by use of a certain Weyl operator U = Wm,n ® 1 with k' = {k + m) mod d 
and /' = (/ + n) mod d, meaning there are no restrictions on translations T in the 
phase space 

■^■■0)^ it?] ■ (3.28) 



J J v. 

There are even more local unitary transformations mapping W onto itself such as 
Un <H) with 

^ d-1 

= ■ (3-29) 

s,t=0 

Under this transformation the Bell states Pk,i become 

Utz ® U^Pk,iUl, ®Ul^ Pi,k ■ (3.30) 
Thus, this is a quarter rotation around the origin in the phase space picture 

7e :(*)_([). (3.31) 

The next transformation we consider is Uy (8> with 

d-1 

Uv = J2w-'^'+'^y^\s){s\ . (3.32) 

This affects the Bell states in the following way 

Uv ® U;Pk,iUl ®U^^ Pk+i,i . (3.33) 
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In the phase space this is a vertical shear 



V : (]) ^ I ^] . (3.34) 



Another realisable mapping is the vertical reflection S 

The corresponding anti-unitary transformation on the Hilbert space TC^^ is the 
tensorial product of local complex conjugation C 

d-l d-l 

|s) — ^^a*|s). (3.36) 

s=0 s=0 

The composite application C ® C of this anti-unitary transformation onto the P^^i 
yields the desired vertical reflection 

C ® CPk,iC^ = P-k,i ■ (3.37) 

The three transformations TZ,V,S together with the translation T are the generating 
elements of an arbitrary phase space transformation of form (see |49] ) 

/) - (; :) (t) + - (; :) • p-^^' 

with det{M) = 1 or det{M) = d — 1. For det{M) = 1 the corresponding transforma- 
tion acting on the Hilbert space Ti,^^ is unitary, for det{M) = d — 1 it is anti-unitary. 
Hence, further transformations as for instance horizontal shear, squeezing, horizon- 
tal or diagonal reflection with varying origins can be decomposed into T,TZ,V and S 
and a representation acting on Ti,^^ can be constructed out of them. Transforma- 
tions that cannot be written in the form (13.381) or that do not obey det{M) = 1 (or 
d — 1) are excluded because they do not possess local ( ant i-) unitary representations. 
A detailed proof of this fact from a group theoretical point of view can be found in 

m- 

Let us briefly point out some major consequences of these equivalences. In most 
cases we study low dimensional sections of the simplex. These slices are mainly 
mixtures of particular Pk/s and the unity 1 which is an equally weighted mixture 
of all Bell states 1 = Ylt~i=o and can be regarded as uncolored noise. The 
transformation rules imply that all one parameter states p = + (yPk,i with 

arbitrary k,l G {0,..,d — 1} but same a have the properties in terms of separa- 
bility and (non-)locality (Identity 1 is mapped onto itself for all unitary transfor- 
mations and the single Bell state P^/ can be translated freely). These are the 
so-called isotropic states. The same applies to two-parameter families of states 
of form p = ^—^2-^1 + aPk^i + l3Pm,n with arbitrary k,l,m,n G {0, ..,d — 1}. We 
translate the flrst point [k, I) to the origin (0, 0) and then bring the second point 

{m — k,n — l) to (0, 1). This can be done by the matrix M = ^ q ^'^^^ 

q{n — l) = 1 = det{M). As a result we have found that such states with different P^^i 
and Pm,n but same a and (3 are equivalent and furthermore they are symmetric in a 
and (3 (i.e. apply the same procedure but interchange the roles of the points). For 
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three-parameter states p = ^—^^f^l + aPk^i + (3Pm,n + lPp,q the situation changes. 
We can proceed as before bringing two points to (0, 0) and (0, 1). If the third point 
was on a hne {(/c, /), (m, n), (p, q) = {a{m — k) + fc, a{n — l) +/)} then it is now on the 
hne with /c = due to hnearity of M, for instance for d=3 this point takes on (0,2). If 
it was not an element of this line it can be brought to (1, 0) with horizontal shear and 
vertical reflection without influencing the points on /c = 0. However, it cannot be 
brought to = without transforming the other two points. This shows that not all 
three parameter states with same a, (3 and 7 necessarily possess the same attributes. 
We do not discuss cases with more than three points in detail but want to stress 
that all complete lines {(/c, /), (m, n), .., {{d — 1)- {m — k) + k, (rf — 1) • (n — /) +/)} are 
equivalent because they can all be mapped onto the hne {(0, 0), (0, 1), .., (0, d—1)} by 
transforming two of the points onto (0, 0) and (0, 1) with the above stated method. 




Figure 3.5: Illustration of all possible complete lines through the point (0,0) for d = 3. Lines can 
be completed by points with k,l ^ {0, ..,d — 1} because of the periodicity of the Weyl operators 
implying Pk,i = Pk+m-d,i+n-d for all n, m e Z 



3.2.3 Geometry of separable and entangled states 

Theoretical strategies for determining the separable and entangled states within W 
are discussed and applied onto sections of the simplex for d = 3. We begin with 
entangled states p that can be revealed by the PPT criterion (see §1.2.21) and provide 
some simplifications for determining definiteness of operators p^-* that are inferred 
from the structure of the simplex W. The explicit form of a projector / is 

^ d-l 

Pk,i = ^J2 \s-l,s){t-l,t\ . (3.39) 

s,t=0 

Therefore the elements of W are 

P=d ^ c,,;w'^^~*)|s-/,s)(t-/,t| (3.40) 

s,t,k,l=0 
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and partial transposition on subsystem A gives 



d-l 

d 

s,t,k,l=0 



P^^^l J2 ^k,iw''^'-'^ \t - I, s) {s -l,t\ . (3.41) 



We substitute the running index I by I — s + t — m which leads to the expression 

d-l 

d 



d-l 

P^^ = -j Ck,s+t-mW^^'-*^ \m - s, s) (m - t, t\ 



s.t,k.'m=0 



d-l 

E 

m=0 
d-l 



d-l 

- ^ Ck,s+t-mw''''''^^ \m- s,s) {m-t,t\ 



s,t,k=0 



= Y,Bm. (3.42) 

m=0 

As we observe, each operator Bm acts on a subspace Tim = C'' spanned by the 
vectors {\m — s) |-s)}s=o,..,d-i- Thus, for determining definiteness of p^-^ it suffices 
to consider the definiteness of the d x d matrices Bm having the entries 

^ d-l 

(^-)s,* = ^E^M+*-m^'^^"*^ ■ (3-43) 

k=0 

Relations between the matrices B^ lead to further simphfications of this problem. 
For example 

^ d-l 

(Bm-2),, ^^Yl Ck,s+t-m+2w'^'-'^ = (5m)s+l,t+l (3-44) 
k=0 

holds for even and odd d and signifies that Bm^2 and B^, are unitarily equivalent. 
Moreover, for odd d all B^. are unitarily equivalent which can be proven by using 
the periodicity (modulo d) 

^ d-l 

k=Q 

d-l 

d 

k=0 

= {Bm)s+{d+l)/2,t+{d+l)/2 ■ 

(3.45) 

Since unitarily equivalent matrices have the same eigenvalues it suffices to compute 
the eigenvalues of a single Bm when the dimension d is odd and two Bm (one with 
odd m and one with even m) when the dimension d is even. 

As we already know, we can exclude states p with non-positive p^^ from being 
separable but positivity p^^ > does not guarantee their separability. We now show 
which states p G W arc separable with certainty. Consider functions kn{x) and ln{x) 
with X G {0, ..,d — l} defining complete lines {(A;„(0), /n(0)), .., {kn{d — l), /„((i — 1))}. 
We state that all states that are equally weighted mixtures of states forming a 
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complete line are separable (in particular these are the outermost separable states 
of W as is shown in [50] ) . 

An = ^$^Pfc„(-)AW e SEP . (3.46) 

x=0 

We have to prove this for one particular complete line only, for instance k{x) = x 
and / = because of their equivalence in terms of local unitaries ( §3.2.2p 



x=0 
d-1 



= ^ E ^^''^'''^\s,s){t,t\ . (3.47) 

x,s,t=0 

Once again we use ^ J2t=o w^^'^~^^ = 5s,t and get 

^ d-1 

-J2\s,s){s,s\ , (3.48) 



d-1 

i 

Pline 



d 

which is obviously separable. All possible convex combinations of all line states A„ 
03.461) form the kernel polytope (KP) 

KP = J p G W I p = J]c„A„ , c„ > , J]c„ = 1 i , (3.49) 

I, n n ) 

a subset of separable states of PPTnVV0. 

At this point the separability of states within the kernel polytope and the non- 
separability of NPPT states are ensured. For the remaining states that are PPT 
but do not lie within the kernel polytope we have to construct optimal entanglement 
witnesses Wopt. if we want to completely clarify the question of separability. As we 
already know for all states p that lie on the boundary of the convex set of separable 
states 9SEP there exists an optimal entanglement witness with Tr{Wopt.p) = 
( §1.2.4p . To determine the boundary states of SEPflW and their witnesses the 
symmetry of the simplex is of great help. Consider a state p which is invariant 
under the symmetry group Q of unitary or anti-unitary operators Vg & Q, meaning 
VgpV~^ = p for all g. Thus, 

TriWp) = Tr{WV-WgpV^Wg) = Tr{VgWV-WgpV;') = Tr{VgWV;'p) 

shows that the symmetries ^ of a state p are reflected in its witnesses. In other words 
we can restrict the search for witnesses on ^-invariant operators VgWV~^ = W (for 
all g). This is a significant restriction on the form of the witness when the regarded 
state p is of high symmetry. In our case where all states belong to the magic simplex 
W all states are invariant under the symmetry group Q = {2Pk i — l}k,i=o,..,d-i and 
the most general form of an ^-invariant W is 

W = Y, ^k,iPk,i fik,i e R . (3.50) 

k,l 



l^When we speak of PPT as a set wc mean PPT={p e Wdf | Tr{p) = 1, p > 0, p^^ > } 
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According to (11.641) W must have non-negative expectation values for all product 
states 



(V^l ® {ri\W\%lj) ® \vi) 



d-l 



2 /tfc,/ (^1 ® {'n\ Wk,i ® 1 |s) O \s) {t\ O {t\ O 1 1^) ® |r/) 

fc,/,s,t=0 
^ d-l 

- ^^k,imwkAs)mm{t\wiM) 

k,l,s,t=0 

' d-l 



.k,l=0 



> 



V 



(3.51) 



Here we have introduced a matrix M<^ that depends on the vector |0) = Ylt=o ivl-^) k) 
which is merely an anti-unitary transformation of \ri) (complex conjugation of the 
expansion coefficients of \ri) in the basis {|s)}). Consequently, if non-negativity 
of W holds for all |0) then it also holds for all and vice versa. Iff W is not 
only an entanglement witness but also optimal, then there exists a product state 
lip') ® W) e so that 



(^'1 ® (VI W ® W) = i (^'1 il/fy 1^') 



. 



(3.52) 



Since is non-negative the eigenvalue for lip') must be zero and if this is the case 
then W is an optimal entanglement witness for p = t]') {tp' , 7]'\ and all incoherent 
superpositions that are compatible with the symmetry of W 



Pg 



c,>0 



(3.53) 



It follows that the boundary of SEPflW is determined by the innermost states of 
the set yy obeying Tr{Wp) = 0, where W is an operator W = Yliki'^k,iPk,i whose 
associated non-negative matrices have at least one vanishing eigenvalue, i.e. 
det{M^) = 0. Finding those innermost states is still a very difficult task even though 
symmetries have narrowed down the search for their optimal witnesses and in many 
cases one must perform a numerical variation of the parameters {Kk,i} and the vector 
10). Analytical solutions can be obtained for states p with further symmetries. For 
instance, if we intend to find the boundary of the one-parameter family of states 



P 



a 



rf2 



1 + aPo,o 



(3.54) 



which is invariant under all phase space transformations except translations, we 
can restrict the search for an optimal witness on = al + 6Po,o having the same 
symmetries. For this witness the associated matrix M<^ is 



M^ = dal + b\ 



(3.55) 



Here the eigenvalues do not depend on |0) and are da + b, da and da so one must 
choose a > and b = —da in order to get an optimal W {det{M^) = 0). The state 
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on the boundary of SEP is then given by 



Tr{Wp) = aTr 



1- a 

[l-dPo,o){—p^l + aPo,o) 



d-1- a{d^ - I) 
d 



= 



(3.56) 



which is achieved for a = It follows that all isotropic states p = + CiPp,q 
have this bound for a due to their equivalence by local unitaries (see §3.2.2p . Their 
optimal entanglement witnesses W = al — daPp^q with a > define the enclosure 
polytope. For an arbitrary state p G W we compute 



Tr{Wp) = aTr 



d-1 



[1 - dPp,g){j2 Ck,lPk,l) 



k,l=0 



d-1 

a{-dcp^q + y^Cfc,;) 

k,l 

a{-dcp^q + 1) . 



(3.57) 



This implies that any p that has at least one component Ck^i > ^ is detected to 
be entangled by one of the isotropic witnesses W = al — daPp^q. Consequently, 
separable states lie within the so-called enclosure polytope (EP) defined by 



EP = <! p G W I - > Ck,i > Vcfe,, 



(3.58) 



We now study two and three dimensional sections of the simplex for d = 3 (qutrits) 
in order to illustrate regions of separable and entangled states that result from these 
concepts. The first family of states we consider is a mixture of two Bell states which 
are all locally unitarily equivalent to p = ^~q~^ l+a-Po,o-l-/3-Pi,o. This means we have 



Co,o 



l-a-(3 



+ a , Cin 



l-a-(3 



+ (3 and all other components are c^^, 



-13 



which 



9 I " ; -i,U 9 

implies a > j3 > and /3 < 1 — a for positivity of p. We obtain the boundary 
of PPT by setting det{Bo) = which yields (2a + 2f3 + l){8a^ + - Ilpa + 2a + 
2/5 — 1) = (details on this and other computations regarding p can be found in 
the Appendix [D]). Due to invariance of p under horizontal reflection we restrict on 
entanglement witnesses of form W = A|l+aPo,o+^-Pi,o + C-P2,o- As it has been shown 
in this section, line states for example pune = |(-Po,i + Pi,i + -^2,1) are separable. 
Hence, Tr{Wpiine) = f > must always be valid and therefore A cannot be negative. 
When we set A = 1 (which only fixes the scaling of the witness) the associated matrix 
becomes = 1 + aWo,o 10) (01 <o + bWi^ |0) (0| + cW2,o (0| Wl^. A 
numerical search for solutions of det{M^) = by varying the parameters a, b, c and 
the vector |0) was done in [51]. It is shown in [52] that there exist optimal witnesses 
Wopt. so that Tr{Wopt.p) = yields the boundaries 4:a'^-5a+A0(3'^ + {17a-14:)(3+l = 
and 4/3^ -5(3 + AOa^ + {17/3 - 14)a + 1 = 0. Graphically this is illustrated in the 
following figures. 
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PPT optimal witnesses 

Figure 3.6: Illustration of the state p = 1 + aPo.o + PPi,o- AH physical states lie within the 

green triangle which represents the border of positivity. The blue lines correspond to the enclosure 
polytope (the outer one) and the kernel polytope (the inner one). PPT states lie within the blue 
ellipse and there is also a small region of bound entanglement (region between the red curve given 
by optimal witnesses and the PPT boundary) 




Figure 3.7: Enlarged illustration of the region of the state p = -!— + aPo.o + f3Pi,o where 
bound entanglement appears (filled red) 
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Next, we consider three-parameter families of states of form p = 1 + aP^^i + 

PPm,n + lPp,q with k,l,m,n,p,q G {0,..,2}. In section §3.2.21 it has been shown 
that for fixed parameters 7 any state of this family is either locally unitarily 
equivalent to the state p = ^"""^"^ 1 + aPo,o + /5i^i,o + 7^2,0 or p = ^"""^"^ 1 + 
aPofi + pPifi + 7-Po,i depending on whether the index pairs {{k,l), {m,n), {p,q)} 
form a line or not. In both cases positivity restricts Q!,/3, 7 to lie within the region 
a > (and all parameter permutations of this term) and j < 1 — a — (3. The 

PPT boundary for states on a line reads (see Appendix [E]) 

(2a + 2/? + 27 + l){8a^ + + 87^ + 2a + 2/5 + 27 - ll/?a - lla7 - 11/^7 - 1) = 
and for states off a line we get (see Appendix [F]) 

- 16a^ - 16p^ - 167^ + 6/?a^ + + + 6p^a + 6/3% + 6^-/^ 

- 12a^ - 12/32 - 127^ + 3/?a + 37a + 3/?7 - 15/?7a + 1 = 0. 

Like in the previous case, the state p = 1 + aPo,o + /3-Pi,o + 7-^2,0 is 

invariant under horizontal reflection and for this reason the search for optimal 
witnesses can once again be restricted on W = |l + clPqa + bPi^ + cP2,o with 

M^ = l + aWo,o 10) (01 + bWi,o 10) (01 + cW2,o |0) (0| Wl^. In [52] a choice 
for the parameters a, b, c in compliance with the constraints for optimal witnesses 
was found that yields 

40a2 + (17/3 + 177 - 14)a + 4/3^ + 7(47 - 5) - /3(197 + 5) + 1 = 
and permutations: (a ^ /3), (a ^ 7), (/? ^ 7) 

for the boundary Tr(Wopt.p) = 0. Unfortunately, for states off a line there exists 
no such solution because obtaining it is much more difficult due to the fact that all 
nine parameters K,k,i of 1^ = Ylk i '^k,iPk,i have to be taken into account because the 
state has fewer symmetries. Regardless of this, for both states the boundaries of 
positivity and PPT are graphically illustrated in the following figures. 
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7 




Figure 3.8: Illustration of the state p = ^-"-^P^i \ _|_ aPo.o + (iPifi + 7-P2,o- Physical states lie 
within the green tetrahedron (positivity) . The boundary of PPT states is a cone (blue) and thus 
all states beyond this surface are entangled. The tip of the cone touches the surface of positivity 
ata = /3 = 7= ^ illustrating the separability of the line state pune 
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PPT 

Figure 3.9: Illustration of the state p = 1+q;Po,o+/?-Pi,o+7-Po.i- Physical states lie within 

the green tetrahedron (positivity). The surface of the PPT state region is a complex geometric 
object (blue surface) 
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3.2.4 Geometry of non-local states 



Our purpose is to determine the set of states p of the magic simplex W that vi- 
olate the CGLMP inequality. For any state of this kind there exists a Bell op- 
erator Bi^ such that 1^ = Tr{pBj^) > 2 and thus a state p for which the in- 
equality maxg^^ Tr{Bj^) < 2 holds is an element of the complementary set. This 
means that the complement and thereof the set itself can be determined by use 
of the optimisation procedure that has been introduced in §2.3.2[ In general, it is 
computationally intensive to determine the boundary with high precision because 
states of W have to be parameterised and varied until a required precision AJ of 
max Id = maxgj^ Tr{Bi^p) = 2 ± AI is reached. However, in the case when for a 
given state yU G W the maximal value of max Ici{p) = max^ j^Tr{Bi^p) is known 
with high accuracy it is not necessary to perform further numerical investigations 
in order to obtain the state on the boundary for the family p = + ap because 
of the relation Tr{BiJ = (see ^ OXTD . 



maxrr(;B/_^p) = max Tr(;B 



Bi, ' Bid ' 



1 — a 

1 + ap 



max Tr{Bi^ap) 

amaxld{p) ■ (3.59) 



Hence, p = -^1+ap lies on the boundary of (non-) locality for a = 2/ max Id{p). For 
our numerical investigations of the three families of states p = ^~Q~^ l-|-ttPo,o+/^-Pi,o, 

P = + aPo,o + f3Pi,o + 7^2,0 and p = ^""/'^ l + «Po,o + /3Pi,o + 7^o,i this 

fact is of great help because out of any numerical obtained value max/3(p) we can 
derive a state on the boundary. More precisely, for a certain choice of a, (3 and 
with resulting value maxJ3(p) it implies that the state with the parameter values 
ab = [2/ max Is (p)] ■ a, Pb = [2/ max/3(p)] ■ /3 and 7b = [2/ max J3(p)] ■ 7 lies on the 
boundary. For instance, the subsequent illustrations of the boundaries of 
(non-)locality were deduced from values max/3(p) of states on the boundaries of 
positivity {a = (and all parameter permutations of this term) and 7 = 1 — a — 

P). For the two-parameter family we calculated max J3(p) for 60 such equally spaced 
points and for each of the two three-parameter families we calculated max 73 (p) for 
920 of them. 



15^=0 for the first family p = ^"""^ 1 + aPo.O + I^Pl,0 



64 



enclosure polytope 
kernel polytope 




POO 



a 



PPT optimal witnesses 



Figure 3.10: Illustration of the state p — ~g~^ 1 + aPo.o + PPi.o- AH physical states he within 
the green triangle which represents the border of positivity. States p in the red filled area violate 
the CGLMP inequahty (max/3(p) > 2) 



The result of the calculations for the state p = ^—^^1 + aPo,o + f3Pi,o is illustrated 
in the above figure. The calculated points on the boundary are connected through 
lines and the area of non-local states is filled red. The boundary of (no n-) locality 
seems to describe a circle for a, /5 > and a line if one of the parameters is negative, 
i.e. q; < or /3 < 0. Suggestions on their specifications are given after the next 
figures illustrating the family p = kz^^S^j _|_ aPo^ + /3Pi,o + tA.o- 
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Figure 3.11: Illustration of the state p = ^— + aPo.o + PPifi + 7-P2,o- States p beyond the 
red/blue shaded surface violate the CGLMP inequality (max/3(/9) > 2). 



7 




non-local states 1 



Figure 3.12: Illustration of the state p — ^ "g^ 1 + q;Po,o+/3J\,o+7-P2,o form a different perspec- 
tive. States p beyond the red/blue shaded surface violate the CGLMP inequality (max/3(p) > 2). 
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The result of the calculations for the state p = + aPo,o + PPi,o + 7-^2,0 

is illustrated in the figures on the previous page. The calculated points on the 
boundary tendentially describe a spherical surface that is intersected by planes in 
the region of negative parameters. Based on symmetries and the facts that the 
boundaries of positivity, PPT and (non-) locality meet at the point a = /5 = 7 = | 
and that the suggested analytical ideal measurements f l2.49p yield the boundary 
parameter |(6v^ — 9) for the isotropic states, we were able to derive the radius 
r = (413v^ — 558) and the center of the sphere a = /5 = 7 = Y|g(— 361 + 186v^). 
These specifications coincide with the numerical data up to the order 10~^ and we 
believe that these discrepancies should decrease for better accuracy and precision 
goals of the numerical optimisation (see §2.3.2p . We suppose that the intersecting 
planes are given by the functions 7 = |(a + /? + 6\/3 — 9) (and permutations a <-> 7, 
/3 ^ 7) because of compliance with the boundaries of the isotropic states and the 
numerical data (also up the order 10^^). It should be noted that the circle and hne 
boundaries in figure 13.101 can easily be obtained on the basis of these specifications 
because it solely illustrates the special case 7 = 0. 

For the remaining family p = + tt-Po,o + f^Pifi + 7-Po,i the geometric form 

of the boundary seems to be more complex and therefore we do not want to make 
any uncovered suggestions on the exact form. To get an impression of this, we have 
illustrated the raw data points in the following figures. Besides the complex shape 
of the boundary, we recognised an interesting peculiarity, namely, in contrast to 
the mixtures of states on a line where there is only one state a = /3 = 7 = |on 
the boundary of positivity 1 — a — /5 — 7 = that does not violate the CGLMP 
inequality, here we have a whole region of local states for 1 — a — /3 — 7 = 0. 
Moreover, in comparison the entire region of non-local states is smaller while the 
region of entangled states is larger. This can be seen as a further example for the 
diverging behaviour of these properties. 
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1 ± POl 




Figure 3.13: Illustration of the state p = ^ "g^ "^ 1 + aPo,o + PPi,o + 7^o,i- R-cd points denote 
states on the boundary of nonlocality (niax/3(p) = 2). 
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Figure 3.14: Illustration of the state p = ^ "g^ ''' 1 + aPo.o + f3Pi,o + lPo,i- Red points denote 
states on the boundary of nonlocality (niax/3(p) — 2). 
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3.3 Geometry of multipEirtite qubit systems (publication) 

In this section we present a paper that was pubhshed on Physical Review A 78 
in collaboration with Beatrix Hiesmayr, Marcus Hubcr, Florian Hipp and Philipp 
Krammer. In this, we discuss a generalisation of the tetrahedron for multipartite 
qubit systems. We investigate separabihty, nonlocality and distillability with regard 
to a multipartite entanglement measure. In the context of this diploma thesis, 
this should be regarded as a further example of a geometric investigation of the 
state space. We emphasise that, because it contains several concepts for the study 
of multipartite systems that have not been introduced. Most of them should be 
self-explanatory or easily comprehensible with the knowledge gained in this thesis. 
Nonetheless, we recommend to read the cited publications. 
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Simplex of bound entangled multipartite qubit states 

B. C. Hiesmayr, F. Hipp, M. Huber, P. Krammer, and Ch. Spengler 
Fakulttit flir Physik, Universitdt Wien, Boltzmanngasse 5, A-1090 Vienna, Austria 
(Received 30 June 2008; published 24 October 2008) 

We constract a simplex for multipartite qubit states of even number n of qubits, which has the same 
geometry concerning separability, mixedness, kind of entanglement, amount of entanglement, and nonlocahty 
as the bipartite qubit states. We derive the entanglement of the class of states which can be described by only 
three real parameters with the help of a multipartite measure for all discrete systems. We prove that the bounds 
on this measure are optimal for the whole class of states and that it reveals that the states possess only «-partite 
entanglement and not, e.g., bipartite entanglement. We then show that this «-partite entanglement can be 
increased by stochastic local operations and classical communication to the purest maximal entangled states. 
However, pure ^-partite entanglement cannot be distilled; consequently all entangled states in the simplex are 
«-partite bound entangled. We study also Bell inequahties and find the same geometry as for bipartite qubits. 
Moreover, we show how the (hidden) nonlocaUty for all «-partite bound entangled states can be revealed. 

DOI: 10.1103/PhysRevA.78.042327 R\CS number(s): 03.67.Mn 



I. INTRODUCTION 

Entanglement is at the heait of the quantum theory. It is 
the source of several new applications as quantum cryptog- 
raphy or a possible quantum computer. In recent years the 
study of higher-dimensional quantum systems and/or multi- 
partite systems has shown that different aspects of the en- 
tanglement feature arise. They may have new applications 
such as multiparty cryptography. 

In this paper we contribute to the classification of en- 
tanglement in a twofold way, i.e., which kind of entangle- 
ment a certain class of multipartite qubit states possesses, 
using the multipartite measure proposed in Ref. [1], and 
whether this kind of entanglement can be distiUed. Our re- 
sults suggest that for multipartite systems one can distinguish 
between different possibilities. 

The class of states we analyze are a generalization of the 
class of states which form the well-known simplex for bipar- 
tite qubits (Sec. n), i.e., aU locally maximally mixed states 
[2,3]. We make an obvious generaUzation and find an analo- 
gous simplex for states composed of an even number of qu- 
bits n, i.e., this class of states shows the same geometry 
concerning positivity, mixedness, separability, and entangle- 
ment (Sec. in). Further, the used multipartite measure [1] 
reveals that the kind of entanglement possessed is only 
n-partite entanglement where n is the number of qubits in- 
volved. The vertex states of the simplex are represented in 
the bipartite case by the weU-known BeU states; for n>2 
they are equivalent to the generalized SmoUn states proposed 
by Refs. [4,6-8]. 

Then we discuss the distillability of the entangled states 
and find states for which the «-partite entanglement can be 
increased by a protocol based only on copy states and sto- 
chastic local operations and classical communications 
(LOCC). We show that the state is not distUlable for any 
subset of parties and hence bound entangled; however, the 
n-partite entanglement can be enhanced to reach the maximal 
possible purity and n-partite entanglement within the class of 
states under investigation, i.e., the vertex states. For a subset 
of these states it has been shown that they allow for quantum 

1050-2947/2008/78(4)7042327(10) 



information concenti^ation (e.g., Refs. [4,5]), so we suggest 
that it tnight stiU be advantageous to enhance the n-partite 
bound entangled states for some appUcations. 

Last but not least, in Sec. VI we address the question as to 
which of the simplex states violate the generalized BeU in- 
equality which was shown to be optimal in this case, and 
draw its geometiical picture. Fig. 4. 

n. THE SIMPLEX EOR BIPARTITE QUBITS 

A single qubit state to exists in a two-dimensional Hilbert 
space, i.e., Ti. = and any state can be decomposed into the 
weU-known PauU matrices cr,-: 

0)= ^(l2 + n,cr,)> 

with the Bloch vector components neR^ and 2^^jn?=|np 
^ 1. For |«p< 1 the state is mixed (corresponding to Tr ii? 
< 1) whereas for |np=l the state is pure (Tr (i?=\). 

The density matrix of two qubits p on C^® C^ is usually 
obtained by calculating its elements in the standard product 
basis, i.e., |00), |01), |10), |ll). Altematively, we can write 
any two-qubit density matrix in a basis of 4 X 4 matrices, the 
tensor products of the identity matrix I2, and the PauU ma- 
trices, 

p = -(I2 ® I2 + afl'i iSi I2 + fo,:l2 ® f,- + CyCr,; (g) (Tj) 

with ai,bj,Cjj e R. The parameters a,- and bj are called local 
parameters as they determine the statistics of the reduced 
matiices, i.e., of Alice's or Bob's system. In order to obtain a 
geometiical picture one considers in the foUowing only 
states where the local parameters are zero (d=b=0), i.e., the 
set of all locaUy maximaUy mixed states, Tr/i(p)=TrB(p) 
= \l2 (see also Refs. [2,3]). 

A state is called separable if and only if it can be written 
in the form 2,p,pf iSipf with p,-3=0, 2p,= l, otherwise it is 
entangled. As the property of separabiUty does not change 
under local unitary ti'ansformation and classical communica- 

©2008 The American Physical Society 
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FIG. 1. (Color online) Here the geometry of the state space of 
even number of qubits is visuahzed. Each state is represent by a 
triple of three real numbers c, Eq. (1). The four black dots at the 
vertices of the cube represent four orthogonal "vertex" states. In the 
case of two qubits these are the foui' maximally entangled BeU 
states i//^,<A* and for higher n they are equal mixtures of 2" I A 
Greenberger-Home-Zeilinger (GHZ) states. The positivity condition 
forms a tetrahedron (red) with the four "vertex" states and the to- 
tally mixed state at the origin (black dot in the middle). AH sepa- 
rable states are represented by points inside and at the surface of the 
octahedron (dashed object). The dashed line represents for n=2 the 
Werner states and for n>2 the generaUzed Smohn states (becoming 
separable when blue changes into green). 

tion, the states under consideration can be written in the form 
[2] 

p=i(l2®l2-FC,-<7,-<8)(7,-), 

where the c, are three real parameters and can be considered 
as a vector c in Euclidean space. Differently stated, for any 
locally maximally mixed state p the action of two arbitrary 
unitary transformations C/j ® C/2 can via the homomoiphism 
of the groups SU(2) and S0(3) be related to unique rotations 
O] igi O2. Thus the correlation matiix Cjj(Tj® <jj can be chosen 
such that the matrix Cy gets diagonal via singular value de- 
composition. Therefore, three real numbers combined to a 
vector c can be taken as representative of the state itself. 

In Fig. 1 we draw the three-dimensional picture, where 
each point c corresponds to a locally maximally mixed state 
p. The origin c=0 corresponds to the totally mixed state, i.e., 
^12®12. The only pure states in the picture are given by 
|c|^=3 and represent the four maximally entangled BeU 
states |i/rt)=(l/v2){|01)±|l0>}, |0±)=(l/v2){|OO>± |ll)}, 
which are located at the vertices of the cube. The planes 
spanned by these four points are equivalent to the positivity 
criterion of the state p. Therefore, all points inside the tetra- 
hedron represent the state space. 



It is well known that density matrices which have at least 
one negative eigenvalue after partial transpose (PT) are en- 
tangled. The inversion of the argument is true only for sys- 
tems with 2® 2 and 2®1> degrees of freedom. The PT cor- 
responds to a reflection, i.e., C2^-C2 with all other 
components unchanged. Thus all points inside and at the 
surface of the octahedron represent all separable states in the 
set. Of course, one can always make the transformation c 
-^-c, and thus one obtains a min'ored tetrahedron, spanned 
by the four other veitices of the cube. Clearly, the intersec- 
tion of these two tetrahedra contain all states which have 
positive eigenvalues after the action of the PT. 

In Ref. [9-11] a generalization to higher-dimensional bi- 
partite states is considered and a so-called magic simplex for 
qudits is obtained. Here the class of all locally maximally 
mixed states have to be reduced in order to obtain this gen- 
eraUzed simplex. Already for bipartite qutrits many new 
symmetries arise and regions of bound entanglement can be 
found (see also Refs. [12-16]). 

We also want to generalize the simplex of bipartite qubits; 
however, in our case we increase the number of qubits. 

m. A SIMPLEX FOR n-PARTITE QUBIT STATES 

Assume we have n qubits. Then a generalization can be 
wiitten as 

P = ^(1 + S c,(7, ® (7,- ® ■ • • ® (7,-) = ^(1 + 2 qaf") . 

(1) 

Obviously, for this generalization we foUow the strategy to 
set the local parameters of aU subsystems j, 
Tri,2,...,;-i,j+i,...,n(p)> to zero, as weU as the parameters shared 
by two parties j,k, Tr,_2,...j_ij+i,...,t_i,«,+i,.„,„(p), zero and so 
on untU n-1 zero. 

Again the state can be represented by a three-dimensional 
vector c. For n=3 the positivity condition pS'O requires' 

kf^i. (2) 

This tums out to be the case for aU odd numbers of qubits 
involved. 

For even numbers of qubits the positivity condition p 
3=0 requires that the vector is within the foUowing four 
planes: 



1 3= 


with 
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+1 










-1 












1-1/ 


1-1/ 



These conditions are exactly the same ones as for the two- 
qubit case n=2, i.e., the four planes above form the magic 
tetrahedron. 



The result is obtained by using a standard computer programm. 

^Note that for n= 4, 8,... the mirrored tetrahedron (c— »-?) is ob- 
tained. Again the result is obtained by using a standard computer 
programm. 
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The purity Tr(p^) gives (I/2")(I + |c|2); thus the states 
with |c|^=3 are the purest states of the class of states under 
investigation and are located in the vertices of the tetrahe- 
dron. Note that with increasing n the percentage of purity 
decreases, i.e., only for n=2 do the vertices present pure 
states. Further analysis of these vertex states foUows later. 

Now we want to investigate if the separabUity condition 
also for n>2 corresponds to the octahedron. The partial 
transpose of one qubit {VT„^ q„|,it) changes the sign in front 
of the erf" matrix, i.e., the y component of the vector c 
changes sign. Therefore the states under investigation are 
entangled by the necessary but not sufficient (one-qubit) 
Peres criterion 



= 3,5,..., |cf=Sl, 



n = 2,4,..., 1- 



(4) 



Taking the partial transpose of two, four, . . .qubits changes 
two, four, . . .times the sign and consequently one obtains the 
positivity criterion (3). Taking the partial transpose of odd 
qubits is equivalent to VTg^ 

For even number of qubits the above Peres criterion im- 
pUes a mirrored tetrahedron, analogously to the bipartite 
case; however, we do not know if the intersection, the octa- 
hedron, contains only separable states. For odd numbers of 
qubits the situation is different and we will not investigate it 
further. 

Now two questions arise: first, are all states represented 
by the octahedron separable and, second, what kind of en- 
tanglement does this class of states possess? 

Let us tackle the second question first. To analyze our 
generalized states p fuilher we use the multipartite entangle- 
ment measure for aU discrete systems introduced by Ref. [1]. 
The main idea is that the information content of any n-partite 
quantum system of arbitrary dimension can be separated in 
the following form: 



/(p) + R(p) + Eip)^ =n 

single property entanglement 



(5) 



where 



I(p)-- 



^) 

single propeity of subsystem s 



(6) 



contains aU locally obtainable information (i.e., obtainable 
information a party can measure on its particle), E(p) con- 
tains all information encoded in entanglement, and R(p) is 
the complementing missing information, due to a classical 
lack of knowledge about the quantum state. The total amount 
of entanglement E(p) can be separated into ra-flip concur- 
rences by rewriting the hnear entropy of all subsystems in an 
operator sum; thus one obtains 



£(p): 



cUp) 



4)(p) +(•••) 



two-flip concurrence three-flip concurrence 



cUp) 



«-flip concurrcnce (7) 

These m-flip concurrences are useful for two reasons: first, 
one can obtain boimds on the operators and thus handle 
mixed states and, second, the authors of Ref. [1] showed (for 
three qubits) that the ra-flip concuixences can be reordered 
such that they give the ra-partite entanglement, which in ad- 
dition coincides with the m-separabrUty [17]. 

Here we extend their result for the states under investiga- 
tion. Due to the high symmetry of the class of states under 
investigation the bounds of the m-partite entanglement can 
be computed and herewith we can reveal the following sub- 
structure of total entanglement E(p), 



E(p)-- 



E{2)(p) 



E0)ip) 



bipartite entanglement tripartite entanglement 

+ ■••+ E^„)(p) 

«-paitite entanglement 

with the subsubstructure 

E(2){p) = E(n)(p) + E(n)(p) + ■■■ + E^^„)(p) + £(23)(p) + 
+ E(2„)(p)+ ■■■ +-E(„-i,„)(p), 



(8) 



E(3)(p) = E(i23)(p)+ •••+£( 



(n-2,/i-r 



,n)(p). 



E{„){p)=E^l2...„){p). 



(9) 



We find that for the states under investigation the only non- 
vanishing entanglement is the n-partite entanglement and it 
becomes (for details, see Sec. IV) 



X maxi 0, ^ max(- 1 + c • n*'',- 1 + ( 



l+c-,7»,-l + c-«^^)) 



(10) 



where X=l except for bipartite qubits, when it is X=2 (the 
reason for this difference is explained later). Hence, we find 
the same condition for being entangled as given by the one- 
qubit Peres criterion. 

Now, if these bounds are exact also for n>2, then aU 
states represented by the octahedron ai'e separable. Indeed, it 
tums out that this is the case. We give the proof of separa- 
bility separately in the Appendix. 

In summary, we have foimd for an even number of qubits 
the same geometry as in the case of bipartite qubits, also 
depicted by Fig. 1. Moreover, we have shown that the mul- 
tipartite entanglement measure proposed by Ref. [1] works 
tightiy as the bountJs are exact and it reveals only n-partite 
entanglement. Let us discuss this result more carefully. 
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For the purest states, |c|^=3, located in the vertices of the 
tetrahedron, the maximal n-partite entanglement results as 
£(„)= 1 except for n=2 when it is -E(„)=2. Thus the amount of 
entanglement for n > 2 is independent of the number of qu- 
bits involved. The reason for the difference can be found in 
the information content of the multipartite system, Eq. (5). 
The maximal entanglement of an n-partite state is n. This is 
the case if and only if the local obtainable information of all 
subsystems is zero and the classical lack of knowledge of the 
quantum state is also zero, i.e., the total state is pure. For 
bipartite qubits, n=2, the vertex states are the Bell states, 
which have maximal entanglement 2 whereas the locally ob- 
tainable infoiTnation S is zero, as well as the lack of classical 
knowledge about the quantum state R=0. 

By construction for n>2 we set the locally obtainable 
information S of all subsystems to zero; however, also all 
possible locally obtainable information shared by two, three, 
. . ., ;i- 1 parties is set to zero; obviously this is not compat- 
ible with being maximally entangled. The information con- 
tent for n > 2 is given by 

n = E„ + R=\+R, (11) 

and consequenfly the lack of classical knowledge is nonzero, 
i.e., R=n-\. Differentiy stated for /j=4, any party has the 
trace state and also any two parties and any three parties 
share the trace state, therefore R='i. 

Remark. The local information S^p) of one subsystem s is 
nothing else than Bohr's quantified complementarity relation 
[18-20], with its well-known physical interpretation in terms 
of predictabihty and visibility (coherence). One can extend 
this concept for two parties sharing a state; then the (bi)local 
information of total multipartite system can be defined in 
similar way and is complemented by the mixedness of the 
shared bipartite system. Again this (bi)local infonnation is 
obtainable only if and only if the state is not the b'ace state. 

Coming back to the simplex geometry we see that the 
closer we get to the origin the more the amount of entangle- 
ment is reduced by increasing the amount of classical uncer- 
tainty R only. 

For bipartite qubits the vertex states |cp=3 ai'e the four 
Bell states. For n qubits we find for |c|^=3 also four unitary 
equivalent states; however, they are no longer pure. For n 
=4 the state is an equally weighted mixture of four |GHZ) 
states. Starting with one GHZ state, e.g., 

|GHZ) = -^{|0000> + |1111)}, (12) 

V2 

one obtains another representation by applying two flips, i.e., 
1®1®(T^® CT^, and then applying on the new GHZ state rep- 



resentation the operator l®cr^®crj®l, and onto that new 
GHZ state representation the operator a^® a^®!®!, giving 
the last GHZ state representation. The other three vertex 
states are obtained by applying only one Pauli matrix. For 
n=6 we have 2* GHZ states where 2^/4 GHZ states equally 
mix for one vertex state. 

Remark. We find the same symmetry for the bipaitite qu- 
bit case, one Bell state is mapped into another by one Pauh 
matrix; however, applying two Pauli matrices maps a Bell 
state onto itself, therefore we have no mixture of different 
maximally entangled states. 

In the next section we give the detailed calculation of the 
measure and in the following section we investigate the 
question whether the entangled states are bound entangled 
and if so in what sense their entanglement is bound. In par- 
ticular we discuss what it means that the substructure re- 
vealed by the measure shows only ;j-partite entanglement. 

IV. DERIVATION OF THE MULTIPARTITE MEASURE 
FOR THE SIMPLEX STATES 

In Ref. [1] a multipartite measure for multidimensional 
systems as a kind of generahzation of Bohr's complementa- 
rity relation was derived. Here, we give exphcitiy the results 
for n=2 and 4 expressed in the farruliar Pauli matrix repre- 
sentation 

It is well known that to compute the concurrence intio- 
duced by Hill and Wootters [21] one has to consider 

p{ay®a,)p*{a,®ay) (13) 

where the complex conjugation is taken in the computational 
basis. The concuo'ence is then given by the formula 

C = max{0,2 max{Xi,X2,^.3,^.4} - (^-i + X2 + ^^3 + ^^4)} 

(14) 

where the X,-'s are the square roots of the eigenvalues of the 
above matrix. To obtain the information content we have to 
multiply this measure by 2. 

The first observation in Ref. [1] is that the linear entropy, 
M(p) = |[l-Tr(p^)] can be rewritten using operators. This 
means, e.g., for any pure four-qubit state, 

1 

\ifi)= 2 aijki\ijkl}, (15) 

ij,k,l=0 

the linear entiopy of one subsystem can be written as 



1 

M\Tr^\,p}{i/^) = M\pi) =2 2 M{(T, ®a,®l® MiMDHMll - [i'Am'iiiklDlrf 

«=0{/,#/;};{,2#,;} 

1 

+ 2 2 Mitr^ ®i®(7^® i){\iiki3i}{iiki3i\ - \i[kg){i[kg\)\if-}f 

{'■|*';};{'3*'3} 
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1 

+ 2 S Mif^x ® 1 ® 1 ® (7j(|;iWi4X;i«;4| - \i[kli^){i[kli:i\)\if}f 

*■■■'=" {i|^^';};{'3'^'3} 
1 

+ 2 S ma ® a-, ® (7, ® I)(,\ki2i3l){ki2i3l\ - \ki:ig){ki2i3l\)\rW 
1 

+ 2 2 Ma ® ® 1 ® crj,)(|A:i2/i4)(yt;2/i4| - |;t/^W;)<fo-y/^|)|i/^)|2 

{i2?^'2!;{'4*'4! 

1 

+ 2 2 Ma ® 1 ® O-^ ® CrJ(|Wl3;4><yt//3/4| - |«i3!4)<yt//3i4|)|lA*)|^ 

'•■•'=0 {i3#/3};{<4'^'4! 

1 

+ 2 2 M(<^x ® (7^ ® (7^ ® l)(|iii2i3/:Xili2l3*:| - |;i;2/3/:>(ii'i2'3^|)|'A'*>l^ 

k {,:,^,-|};{,-2^t,-^};{,3^t,-5} 

1 

+ 2 2 M(<^x ® ® 1 ® ai)(|iii2fa4><'l'2^i4l - l'l'2^'4><'l'2^i4l)l^'^>P 

{,:|^,-;};{,-2^t4};{,4#,;} 

1 

+ 2 2 M(<^x ® 1 ® (7^ ® (7.,)(|/i/:l3l4>(iii:!3!4l - |/;fa3/4)(/;fcj3/4|)| ^4*)^ 

« {/:i?'<-;}*3*--^};{-4?'-|} 
1 

+ 2 2 Ma ®<Tx®(T^® a,){\ki2hU){ki2i3U\ - l^'2'3'4)<^'2'3'4l)l^^)P 

*=0 {/2*/2}*3*'3};{'4*'4> 

+ 2 KiAKo-x ® (7, ® (7, (g) (7j(|/i/2(3/4>(;ii2;3i4l - l'Ii2'3'4)('I'2'3'4l)l'A*)l^' (16) 

{'1 * ' I l^'l*^ 'ii'fc '3}'{'4'^ 

where, e.g., {;i}=i^{;J}, {12} =^{'2} means that the set of indices are not the same, i.e., the sum is taken over 

{;1,;2;;;,J2} = {0,1;0,0},{0,0;0,1},{0,1;1,0},{0,0;1,1},{1,1;0,0},{1,0;0,1},{1,1;0,0},{1,0;0,1}, 

{0,0;1,0},{1,0;0,0},{0,0;1,1},{I,0;0,I},{0,1;1,0},{1,1;0,0},{0,1;1,1},{1,1;0,1}. (17) 



Likewise the linear entropies for the other subsystem can be 
derived, i.e., separated into terms where the flip operator a;^ 
is appUed two, three, or four times. It is well known that for 
pure states the sum over the entropies of all reduced density 
matrices is an entanglement measure; therefore using the Un- 
ear entropy we get the following entanglement measure: 
4 4 
Em - 2 M\p,) = 2 [Cn-A)]', (18) 

s=l m=2 

where (C™)^ is the sum of all terms of all reduced matrices 
that contain m-flip operators. These quantities were called 
(squared) ra-concurrences, because they play a similar role as 
the Wootters concun'ence. 

For mixed states p the infimum of all possible decompo- 
sitions is an entanglement measure 

E{p) = inf 2 PiEm). (19) 

p/.k',) Pi.],/,,) 

The problem of the whole entanglement theory is that this 
infimum can in general not be calculated. Now we bring the 



operator representation of the hneai' entropy into the game, 
because for operators upper bounds can be obtained. 

Let us start with the calculation of the four-flip concur- 
rence C*"*', which is the sum of all terms containing four-flips 
of the entropies of all reduced matrices, i.e., 

[C^'Kp)f= inf 2 pLC^'K^Pi)?. (20) 

pM /),-,k,> 

As shown in Ref. [1] one can derive bounds on the above 
expression for any m-flip concurrence by defining, in an 
analogous way to the HiU and Wootters flip density matrix 
[21], the ra-flip density matrix 

P:'=oM,MQ\-\{iMQ\)p* 

><OMMi,,}\-\{i'MQ\) (21) 

and calculating the X,'„'s which are the squared roots of the 
eigenvalues of p"p. The bound B*'"' of the m-flip concur- 
rence C5"'' is then given by 
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B"'{p) - ( 2 max[0,2 max({\;,}) - 2 



n-partite Entanglement 



(22) 

From Eq. (16) we see that for the four-ffip concurrence of 
subsystem pi four different operators occur, thus we have in 
total 16 diiierent operators Usted in Appendix B. 

Inserting our class of states we find that for each operator 
C the eigenvalues are the same, i.e., one obtains eight zeros 
and the remaining four eigenvalues are exactly equivalent to 
the Peres criterion Eq. (4). 

The same procedure has to be appUed to calculate the 
three-flip and the two-flip concurrence. As can be seen from 
Eq. (16) here the unity and cr- matrices are involved which 
lead to no contribution for the states under investigation. 
Remember that they are mixtures of the vertex states, which 
are equal mixtures of such GHZ states that differ by two 
flips. Therefore the total entanglement is given by the C''*' 
concun'ence only and is a four-partite entanglement. For n 
= 6,8,... the scenario is the same, because of the same un- 
derlying symmetry. 

In Appendix A we show that all states not detected by the 
measure are separable; thus the bounds are optimal and 
therefore the measure detects all bound entangled states. 

V. ARE THE ENTANGLED STATES BOUND 
ENTANGLED? 

In Refs. [4,6-8] the special states c=Ci=-C2=Ct, for n 
>2, which were named generaUzed Smolin states (for n=2 
these states are the Wemer states), are investigated and they 
show that for ls=c> j these states are bound entangled. In 
particular, the authors argued that these states are bound en- 
tangled, because the states are separable against bipartite 
symmetric cuts hke 12|34. .., 14|23....... and therefore no 

Bell state between any two subsystem can be distilled. This 
is obviously also the case for the whole class of states imder 
investigation. 

As the considered measure of entanglement revealed only 
n-partite entanglement and, e.g., not, ra-partite entanglement 
{m<n), it may not seem directly obvious that Bell states 
(bipartite entanglement) cannot be distilled, because the class 
of states does not possess any bipaitite entanglement. Thus 
the question could be refined to ask whether n-partite pure 
entanglement can be distilled. 

For the ^-partite class of states imder investigation we 
consider a similar distillation protocol as the recurrence pro- 
tocol by Bennett et al. [22]. For that we generaUze it such 
that each party gets a copy onto which a unitary bilateral XOR 
operation is performed and afterward a measurement in, say, 
the z direction, is peiformed. Only states are kept where all 
parties found their copy qubit in say, the, up direction. This 
protocol favors, as do, all protocols one state; in our case for 
n=2 it is the <I>+ state and for n> 2 its equivalents. 

In detail it goes as follows. We consider one state and its 
copy 
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-{p"+c,(7r} 



(23) 



and all parties get a copy state. Therefore, we reorder the 
state by a unitary transformation such that the first and sec- 
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FIG. 2. (Color online) «-partite entanglement of the Wemer 
states n-2 (here the y axis has to be multiplied by 2) or the gener- 
alized Smolin states n > 2 before and after the application of the 
introduced protocol (upper dashed green curve). Note that the ver- 
tex states are mapped onto themselves by the given protocol. 

ond terms in the tensor product belongs to Alice and the third 
and fourth terms to Bob, and so on: 



{2") 

+ C/c/cr, ( 

Now each party performs on its two subsystems a unitary 
XOR operation 



[(1 CD I)®" + c,<l «> (T-,-)®" + q(ai ® 1)®" 

Tjf"]. (24) 



1 




\o 





1 

0/ 



(25) 



and then projects on the copy subsystem with P=2(l+cr,). 
This gives again a state in the class of states imder investi- 
gation, i.e., one finds 



l + cl 
2c- 



(26) 



Comparing with the separabUity condition and with the posi- 
tivity condition, one verifies that only separable states are 
mapped into separable states. 

Let us consider the Wemer states and the generaUzed 
SmoUn states (c=c^=c^.=c^), for which we derive that the 
n-partite entanglement is always increased after the above 
protocol (see Fig. 2). For -l/\3=Sc=s 1/3 the measure be- 
fore and after the protocol is zero and for c= 1 the state is 
mapped onto itself. For l/3<c< 1 the entanglement of the 
distilled state is increased compared to the input state. In Fig. 
3(a) we give the three-dimensional picture of how the initial 
state c=0.5 moves after each step toward the vertex state. 
Note that the states are no longer in the set of the generaUzed 
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Mixedness 



FIG. 3. (Color online) (a) Shows the final states after each step of the introduced protocol of an inifial Werner or Smohn state c=0.5, 
where each (green) point represents the obtained state after one step of the protocol, (b) Shows the mixedness [2"/(2"- 1)][1 -Tr(p^)], versus 
«-partite entanglement diagram (for n=2 the y axis has to be multiplied by 2), where the (blue) curve corresponds to the Werner or Smohn 
state whereas the (red) curve is the state connecting two vertices. All states of the simplex have their mixedness-entanglement ratio between 
these two curves. The middle (dashed, green) curve corresponds to the final states of a distilled Werner or Smohn state, and the (green) points 
represent the final states after each step of an initial Werner or Smolin state c=0.5. 



Smolin sets, another advantage of the considered set of states 
as no random bilateral rotation to regain the rotational sym- 
metry is needed. In Fig. 3(b) we show the mixedness- 
entanglement relation of this example. Note that all states of 
the simplex are within the two curves and the middle curve is 
the result for the generaUzed Smolin state after one step of 
the protocol. 

Remark. Not aU states of the simplex are mapped into 
more entangled states by this protocol. For example, the mix- 
ture of two vertex states [c^=(0,0,c) with c=^ 1] is left in- 
variant. 

In summary, we have found a protocol that increases the 
amount of entanglement with local operations and classical 
communication only and the final states are always within 
the class of states. Only for n=2 the final state is pure and 
maximal entangled and therefore the above protocol is a dis- 
tillation protocol, i.e., pure maximally entangled states can 
be obtained. However, for n > 2 the final state is no longer 
pure, but has the maximal n-partite entanglement of the class 
of states under investigation. 

Thus the next logical step is to search for a distillation 
protocol which distills the veitex states into pure maximally 
entangled states, i.e., GHZ states. However, this is not pos- 
sible for the following reasons. In general, any equally 
weighted mixture of two maximally entangled states cannot 
be distilled by mainly two observations. As for all maximally 
entangled states p,- obviously the entanglement can be re- 
duced only by any completely positive map A:p/i-^p,', i.e., 
ii(p')=S£'(p,)VA. And as the entanglement E{p) is convex, 
i.e., E(p'i)+E(pj)'i2E(p'i), we conclude that at least one p,- 
must be mapped unitary onto itself or another maximally 
entangled state. Because aU maximally entangled states are 
equivalent by local unitaries, such a map consequendy maps 
also the other maximally entangled state of the mixture into 
a (different) maximally entangled state. Hence, for no 
equally mixture of maximally entangled states a maximally 
entangled state can be distilled. Note that in the case of bi- 
partite qubits this is trivially true, because any equally mix- 
ture of BeU states is separable, however, for multipartite 
states this is not necessarily the case (e.g., our vertex states). 

Thus we find that we can increase the amount of the 
^-partite entanglement until the vertex state, but not further- 



more, and therefore aU entangled states are bound entangled, 
i.e., no pure «-paitite entanglement can be distilled among 
any subset of parties using stochastic LOCO. The common 
definition of distillation is that no piu'e maximally entangle- 
ment among any subset of parties using LOCC can be ob- 
tained; see, e.g., [23,24]. A different way to prove that the 
entangled states are bound is given in Ref. [25], where it is 
shown that if no singlets can be distilled also no GHZ state 
can be obtained. Therefore for the class of states under in- 
vestigation we also can not distill any bipartite entanglement. 



VI. THE GEOMETRY OF THE STATES VIOLATING THE 
CHSH-BELL INEQUALITY 

Analog to the bipartite qubit state one can derive a 
(Qauser, Home, Shimony, Holt) CHSH-BeU type inequahty 
for n qubit states [26]. Here n- 1 patties measure their qubit 
in the direction a or a and the ;jth party in the direction b or 
b\ then one obtains the following BeU inequahty: 

Tr(SBeU-CHSHP)=S2 (27) 

with 

SgeU-CHSH = ® 'J'? ® ••• ® ®{b + b')a 
n-l 

+ d'a®d'(T » ... ®a'a ®{b-b')a 
\ ^ / 

n-l (28) 

where a,a,b,b' are real unit vectors and the value 2 is the 
upper bound on any local reahstic theory. 

It is known that for n=2 the maximal violation by quan- 
tum mechanics can simply be derived by the state p itself 
[27]. A matrix p violates the BeU-CHSH inequahty if and 
only if M(p}>l, where M(p) is the sum of the two largest 
eigenvalues of the Hermitian matrix OC with (C)y=Tr((7,- 
® Cjp). A generalization for n qubits is simple, because the 
matiTx C is diagonal for the states under investigation, and 
thus the same proof works. 
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FIG. 4. (Color online) The three cylinders show the saturation of 
the Bell inequahty. All states outside these cylinders violate the Bell 
inequahty._The vertex states violate the Bell inequality maximally, 
i.e., by 2\2. 

In our case M(p) is simply the sum of the two largest 
squared vector components. In particular, if ci and C2 are 
greater than we obtain the Bell inequality 

c? + ci=sl. (29) 

This gives a simple geometric interpretation of aU states vio- 
lating the Bell inequahty. AU possible saturated BeU in- 
equahties give three different cyhnders in the picture repre- 
senting the state space (see Fig. 4). AU states outside of these 
three cyhnders violate the BeU inequahty. 

Furthermore, this result shows that an entangled state not 
violating the Bell inequahty (27) can be transformed via the 
introduced protocol into a state violating the BeU inequaUty, 
leading to the conclusion that aU entangled states of the pic- 
ture have nonlocal features. Moreover, in agreement with 
Ref. [28], the possibihty to construct reahstic local models or 
not is no ciiteiion for being bound entangled or not. 

Let us also remark that^Wemer states (n=2) violate the 
BeU inequaUty for c> 1/V2 whereas successful teleportation 
requires only c> 1/2. 

Vn. SUMMARY AND DISCUSSION 

We generaUzed the magic simplex for locaUy maximally 
mixed bipartite qubit states such that we add even numbers n 
of qubits and set aU partial traces equal to the maximally 
mixed states, i.e., no local information obtainable by any 
subset of parties is available. This class of states can be 
described by three real numbers, which enables us to draw a 
three-dimensional picture. Interestingly, we find the same ge- 
ometry concerning separabihty, nuxedness, kind of entangle- 
ment, amount of entanglement and nonlocahty for aU even 
numbers of qubits (see also Figs. 1 and 4). 

For n>2 the purest states, located in the vertices of the 
simplex, are not pure except in the case of bipartite qubits 



(n=2). We show how to derive a recently proposed measure 
for all discrete multipartite systems [1] in this case. For 
mixed states only bounds exist, however, we show that they 
are for the class of states optimal by proving that aU states 
not detected by the measure are separable. 

The measure reveals that these states possess only 
n-paitite entanglement and no other kind of entanglement, 
e.g., bipartite entanglement. The information content of the 
states can be quantified by the generalized Bohr's comple- 
mentarity relation for n>2 

n = S+E„ + R=l+R, (30) 

where R is lack of classical knowledge and 5=0 the local 
information obtainable by any party. 

Then we investigated the question whether the n-partite 
entanglement can be distiUed. We find a protocol using only 
local operation and classical communication which increases 
the «-paitite entanglement to the maximal entanglement of 
the class of states under investigation. These states are the 
vertex states of the simplex; for n=2 they are the BeU states 
and for n>2 they are equal mixtures of such GHZ states 
which are obtained by applying only two flips, cr,.. 

For bipartite qubits n=2 this protocol is a distiUation pro- 
tocol, i.e., pure maximally entangled states are obtained. For 
n > 2 the vertex states are not pure, therefore we search for a 
distiUation protocol that leaves the class of states under in- 
vestigation to obtaui a pure n-partite maximaUy entangled 
state, i.e., the GHZ states. Indeed, we argue that such a pro- 
tocol cannot be found; more precisely, any equal tnixture of 
GHZ states cannot be distUled. Thus for the class of states 
under investigation aU entangled states are bound entangled 
and herewith we found a simplex where aU states are either 
separable or bound entangled. 

In detaU, we show how an initial state moves after each 
step of the protocol increasing the entanglement in the sim- 
plex (see Fig. 2). Moreover, we find that the states violating 
the CHSH-BeU-like inequality, which was shown to be opti- 
mal in this case, have for aU even numbers of qubits the 
same geometry (see Fig. 4). These two results taken together 
mean that one can enhance the n-partite bound entanglement 
by using only LOCC untU the BeU inequaUty is violated. 
Therefore, for aU -partite bound entangled states its (hid- 
den) nonlocahty is revealed and in agreement with Ref. [28] 
the possibihty for a local reaUstic theory to be constructed is 
not a criterion for disfiUabUity, and likewise whether its en- 
tanglement can be increased by LOCC is also no criterion. 

Our results suggest also that one can distinguish between 
bound states for which a certain entanglement measure can- 
not be increased by LOCC (in our case the vertex states) and 
states for which the entanglement can be increased by 
LOCC, which may be denoted as "quasibound" entangled 
states (aU bound entangled states of the class except the ver- 
tex states). The introduced (distiUation) protocol distUls 
maximally entangled states within the set of states which are, 
however, not pure, but the purest of the set of states. 

Last but not least we want to remark that a subset of the 
class of states has been considered in the hterature, e.g., 
[4,6-8], the so-caUed SmoUn states. For which it was shown 
that no BeU states may be distUled. The theorem in Ref. [25] 
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states that if and only if bipartite entanglement can be dis- 
tilled then also GHZ states — in our tenninology /j-partite 
entanglement — can be distilled. 

In summaiy, we have shown in this paper expUcitly that 
the multipartite measure proposed by [1] detects all bound 
entanglement in the class of states and that the states do not 
possess bipartite entanglement and how the ;j-paitite en- 
tanglement can be increased to a certain value. 

These results not only help to reveal the mysteries of 
bound entanglement by refining the kind of entanglement, 
but they may also help to construct quantum communication 
scenaiios where bound entangled states actually help to per- 
form a ceitain process [29]. This is clearly important, when 
one has future appUcations in mind, e.g., a multipartite cryp- 
tography scenario. 
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APPENDIX A: PROOF THAT ALL STATES 
REPRESENTED BY THE OCTAHEDRON ARE 
SEPARABLE 

To prove that all states represented by the octahedron are 
separable, we show that this is the case for the following 
points in the octahedron 

















?) 











As any convex combination of separable states have to be 
also separable, we have finalized the proof. We start with n 
= 2 and show how this construction generalizes for n 
=4,6,.... 

Suppose Ahce prepares her qubits in the following two 
states: 

<.=^a2±^<7,), (A2) 

where r, is a Bloch vector pointing in the ; direction and is 
given by any number in [-1 , 1]. Bob does prepares his qubits 
in the very same way. If AUces chooses the positive ; axis 
Bob does the same, and if Ahce chooses the negative sign. 
Bob does the same, thus they share the following separable 
state if the preparation is done randomly with the same prob- 
abiUty: 

(A3) 

These states represent three vertices of the octahedron, thus 
the proof is finahzed for n=2. 



Exphcitiy, we find that for the generahzed Smolin state 
(ci=C2=C3=c), the following state is derivesd: 

ft. = 2 3 Pt' = 4 (l4 + 2 « <r}j , (A4) 

therefore as rf-rf^l-l,!] the generahzed Smolin state is 
separable for p e [-^ , ^]. 

For n=4 we remark that with the combination 

= ® + ® = - ■ i^fTi ® a-,) 

(A5) 

one obtains the minus sign, and for the very same construc- 
tion Ahce, Bob, Charhe, and Daisy obtain the following 
separable states: 

AB ^ AB ^ CD , ^ AB ^ CD 

= ^(l4 + rf-rf-rf- cr,- (g) cr,- (g) cr, ® cr,) . {A6) 

As the combination +-, -+ gives again the minus sign, 
this proof generahzes for any even n. 

APPENDIX B: ALL FOUR-FLIP OPERATORS FOR n=4 

For convenience of the reader we Ust all four-flip opera- 
tors in the Pauh matrix representation: 

O' = -{cr,, ® (7,, ® (7y iSl ay - ay ® ay <S) CTj ® (Jj - ay ® CTj 

ay ® <Tj - ay ® cr, cr, ay}, 

= -{cTy ® ay ® ay® ay- ay ® ay® cr, ® cr, + cr,, ® cr, 
® cr,, ® cr^ + cr,, ® cr, ® cr, ® cr,,}, 
1 

O' = -{ay ® ay ® ay ® ay + ay ® ay ® cr, ® cr, - cr,, ® a^ 
® cr,, ® cr^. + cr,, ® cr^ ® cr^ ® cr,,}, 

O* = -{cr,, ® ay ® ay ® ay + ay ® ay ® cr, ® cr, + cr,, ® cr, 
® cr,, ®crj-cr„®cr.,® cr,® cr,,}. (Bl) 

= -{ay ® cr,, ® cr,, ® cr,, - cr^ ® cr,, ® cr, ® cr^ - cr^ ® cr^ 
® cr,, ® cr^ - ay ® ay ®a^® crj, 

Cfi = -{ay ® cr,, ® cr,, ® cr,, - cr^. ® cr,, ® cr^ ® cr,, + cr, ® cr,, 
® cr,, ® cr^ + cr,, ® cr,, ® o\. ® o\.}. 
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O' = ^{fy ® o-y® a-y® ay + (7^0 (Ty® a-^® a-y- a-^® a-y O'^ = ^{cy ^ Cy'^ o'y'S) a-y + cr^® cr^® a-y® a-y + cr^® (Ty 
® ay® ay® ay® a^® a^, ® ay® a^- ay® a^^® ay® a^}, (B3) 

= -{ay ® ay® ay^ ay + a.^® ay® a^® ay + a.^® ay O^^ = -{ay ® ay® ay® ay- a^® a^® ay® ay- a^® ay 
® ay ® aj^ - ay ® ay ® ajc ® a^. (B2) ® a-^® ay - ay® a-^® a-^® oy}, 

= -{o\, ® ay® ay® ay - a,^® a^® ay® ay - a,^® ay O''' = -{cTj, ® ay® ay® ay- a^ ® a,. ® ay® ay + a^® ay 
®ay® a^^-ay ®a^®ay® aj, ® a^^® ay + ay® a^^® a^^® ay}, 

O^" = ~{ay ® ay® ay® ay- a^® a^® ay® ay + a^® ay O^^ = ~{ay ® ay® ay® ay + a^® a^® ay® ay- a^® ay 
® ay® a^ + ay® a^^® ay® a^}, ® a^^® ay + ay® a^^® a^^® ay}. 



O^^ = "{cy ® ay® ay® ay+ a^® a^® ay® ay — a^® ay O'* = ~{o-j, ® ay® ay® ay + a^® a^® ay® ay + a^® a^ 
ay® aj^ + ay® a,^® ay® a^^, ® a,^® ay- ay® a,^® a,^® a^, (B4) 
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Summary and outlook 



It is one of the most seminal discoveries that quantum physics contradicts local re- 
alism. This characteristic can be ascribed to quantum entanglement and manifests 
itself through the violation of a Bell inequality. Closer considerations within recent 
years have raised the question of whether there is a discrepancy between entangle- 
ment and nonlocality. The clarification of this question is of great importance for 
our understanding of the theory, and also might have consequences for future appli- 
cations. The aim of this diploma thesis was to give an overview on the latest state 
of knowledge in this field, and to confront entanglement with the violation of the 
CGLMP inequahty for a certain set of states that is called the magic simplex. 

We began our study in Chapter 1 with an introduction to the mathematical frame- 
work including important operator bases such as the Pauli operator basis for qubits 
and its generalisations for d-dimensional systems, that is the Gell-Mann and the 
Wcyl operator basis. Afterwards, wc discussed several issues concerning entangle- 
ment. Wc showed that it is very demanding to establish whether a mixed bipartite 
qudit state with ci > 2 is separable or entangled, because of the fact that all known 
practical criteria (e.g. PPT and matrix realignment) are only necessary but not 
sufficient for separability. Wc reviewed the requirements for entanglement measures 
and presented possible candidates based on distance and convex roof. In addition, 
we provided the proof that entangled PPT states cannot be distilled via "local op- 
erations and classical communication" . 

In Chapter 2 we made clear the contradiction of quantum physics with local realistic 
theories. In particular, wc explicitly derived that any local realistic description of 
a bipartite d-dimensional system obeys the CGLMP inequality, whereas quantum 
physics predicts a violation for the maximally entangled state. We investigated the 
problem of determining if a given state p is able to violate the CGLMP inequality 
and therefore non-local. This is a high-dimensional nonlinear optimisation problem 
where all possible measurement settings of both parties have to be taken into ac- 
count. An analytical solution has been found only for the case d — 2. In order to 
study the nonlocality of higher dimensional systems, we have developed a numerical 
optimisation algorithm which utilises the generalised Euler angle parameterisation 
of SU{N) and the Nelder-Mead method. The advantages of our method are its 
robustness against local maxima and its adaptability to other Bell inequalities. The 
chapter concluded with an overview of recent approaches for clarifying the possible 
discrepancy between nonlocality and entanglement, such as a discTission on hidden- 
nonlocality, tight Bell inequalities, the maximal violation of the CGLMP inequality 
through non-maximally entangled states and the non-local machine. 

In Chapter 3 we studied the state space of bipartite qubits with the aim of find- 
ing simple representatives of locally unitarily equivalent states. We showed that 
all locally maximally mixed states can be represented by elements of a tetrahedron 
spanned by the Bell states. We discussed the difficulties with regard to finding an 
extension for qudits. We introduced a possible generalisation in form of the magic 
simplex, which is a mixture of maximally entangled two-qudit states. Afterwards 
we considered the problem of separability for this set of states, and showed how it 
can be simplified by exploiting symmetries. Considering two- and three-parameter 
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families of qutrit states {d = 3), we used these concepts to determine PPT states and 
to construct optimal entanglement witnesses. In order to reveal the non-local states 
of this families, we applied our numerical optimisation algorithm. As expected, we 
found that there is a large region of entangled states that do not violate the CGLMP 
inequality. We made a supposition about the exact form of the boundary of CGLMP 
violation for the three-parameter family of states on a line, which coincides with the 
numerical data up to the order 10^^. Comparing two three-parameter families, we 
revealed that for the mixtures off a line, the region of nonlocality is smaller while 
the region of entanglement is larger. This has demonstrated that entanglement and 
nonlocality do not behave conformably in any case. In addition we presented a pub- 
lication of our group, where we investigated a generalisation of the tetrahedron for 
multipartite qubit systems. Investigating separability, we realised that all entangled 
states in this tetrahedron are bound entangled. We showed that there is a relation 
between distillability and the available type of entanglement. In the multipartite 
case, entanglement can be shared in many different ways. In our publication we 
argued that bipartite Bell states cannot be distilled from the occuring entangled 
states, due to the fact that they only possess n-partite (n 7^ 2) entanglement. 

We conclude with some remarks regarding future research. As we have seen in this 
thesis, there is an abimdancc of open problems that make further considerations 
desirable. In order to solve them, progress is needed in the theory of nonlocality 
and entanglement. This means that Bell inequalities have to be improved, or alter- 
natively it has to be shown that this is impossible. In addition, it is necessary to find 
more advanced techniques to solve the separability problem. Without developments 
in this field high- dimensional and/or multipartite systems are almost impossible 
to study. Once this has been achieved, we might be able to fully understand the 
relation between nonlocality and entanglement. 
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A MATHEMATIC A notebook: Partial transposition of mul- 
tipartite density matrices 



ln[1]:= basis [diinension_] : = 

Do[basevector[n] = Transpose [{UnitVector [dimension, n]}], (n, 1, dimension}] 

ln[2]:= register [dimension_, systems_] : = 

Do [eniimeration[n, a_] = Floor [Mod[ (a - 1) / dimension'^ (n - 1) , systems]] +1, {n, 1, systems}] 

ln[3]:= registersHap[dimension_, systems_, s\ibsystem_] : = 

Do[If[n systems + 1 - svibsystem, enumerationswap[n, a_, b_] = 

Floor [Mod[ (a - 1) / dimension"^ (n - 1) , dimension] ] + 1, eniimerationswap[n, a_, b_] = 
Floor [Mod [ (b - 1) / dimension'^ (n - 1) , dimension]] +1], {n, 1, systems}] 

ln[4]:= TensorProducts [systems_, a_] : = 

{V = {1}; Do[v = KroneckerProduct[basevector[en\raieration[n, a]], v] , {n, 1, systems}]; v} 

ln[5]:= TensorProductsswap[systems_, a_, b_] := {v= {1}; 

Do[v = KroneckerProduct[basevector[eniimerationsHap[n, a, b] ] , v] , {n, 1, systems}]; v} 

ln[6]:= pt[matrix_, systems_, subsystem_] : = 

{dimension = (Dimensions [matrix] [ [1] ] ) (1 / systems) ; basis [dimension] ; 

register [dimension, systems]; registerswap [dimension, systems, subsystem]; 
Sum[ TensorProductsswap[ systems, xl, x2] [ [1] ] . 



ConjugateTranspose[TensorProductsswap[ systems, x2, xl] [[1]]] *matrix[[xl, x2]], 
{xl, dimension'^ systems} , {x2, dimension^ systems} ]}[ [1] ] 



ln[7]:= (* partialtranspose for multipartite systems of same dimension *) 
ln[8]:= (* enter: pt [' matrix ',' systems ',' subsystem ' ] *) 

ln[9]:= (* matrix: self-explanatory, systems: number of systems (Integer), 
subsystem: target system for transposition (Integer) *) 

ln[10]:= (*Example*) 



/I 2 8 2 



8 3 2 4 




ln[11]:= pt 



9 6 12 



3 5 i 1 I 



n 8 8 2 ^ 



2 3 2 4 



Out[11]= 



9 3 14 
,6 5 2 7, 
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B MATHEMATICA notebook: Optimisation of Bi^ 



ln[l]:= (* The Gell-Mann matrices are defined *) 








-i 


1 




■ 1 













-i 1 




1 








ln[2]:= A2 = 


i 








; A3 = 





-1 





; A5 = 








; A8 = 1/Sqrt[3] 





1 







. 










, 





, 




, i 







, 





-2 , 



ln[3]:= (* The generalized Euler angles for the four observables A1,A2,B1,B2 are introduced *) 

ln[4]:= QAl = Array [oAl, 8]; Element [QAl, Reals] 

Out[4]= (aAlCl) I q;A1(2) | aAl(3) | aAl(4) | q;A1(5) | q;A1(6) | q;A1(7) | q;A1(8)) e R 

ln[5]:= QA2 = Array [aA2, 8]; Element [QA2, Reals] 

Out[5]= (q;A2(1) I aA2{2) | aA2(3) | aA2(4) | q;A2(5) | q;A2(6) | q;A2(7) | q;A2(8)) e R 

ln[6]:= nBl = Array [oBl, 8]; Element [OBI, Reals] 

Out[6]= (aBUl) I q:B1(2) | q:B1(3) | q:B1(4) | aBl(5) | q;B1(6) | q;B1(7) | q;B1(8)) e R 

ln[7]:= QB2 = Array [aB2, 8]; Element [QB2, Reals] 

Out[7]= (a;B2(l) | q;B2(2) | q:B2(3) | q;B2(4) | aB2(5) | q;B2(6) | q;B2(7) | q;B2(8)) g R 

ln[8]:= AV = Flatten [{QAl, QA2, QBl, £3B2}]; 

ln[9]:= Do[NB[i] = {AV[ [i] ] , 0, 2*Pi}, {i, 1, 32}]; domain = Array [NB, 32]; 

ln[10]:= (* The Euler angle parametrizations of UAl , UA2 , UBl , UB2 are defined *) 

ln[11]:= UAl :=MatrixExp[l*A3*QAl[[l]]] .MatrixExp[l*A2*QAl[[2]]] . 

MatrixExp[i*A3*nAl[ [3] ] ] .MatrixExp [i * A5 * QAl [ [4] ] ] .MatrixExp[i * A3 * QAl [ [5] ] ] . 
MatrixExp[i*A2*QAl[ [6]] ] .MatrixExp [i * A3 * QAl [ [7] ] ] . MatrixExp [i * AS * QAl [ [8] ] ] 

In[l2]:= DA2 : = MatrixExp[i * A3 * QA2 [ [1] ] ] . MatrixExp [ i * A2 * QA2 [ [2] ] ] . 

MatrixExp [i * A3 *QA2 [ [3] ] ] .MatrixExp [i * AS * QA2 [ [4] ] ] . MatrixExp [i * A3 * QA2 [ [5] ] ] . 
MatrixExp [i * A2 *QA2[ [6] ] ] .MatrixExp [i * A3 * QA2 [ [7] ] ] . MatrixExp [i * A8 * QA2 [ [8] ] ] 

ln[13]:= UBl := MatrixExp[l*A3*QBl[ [1] ] ] .MatrixExp[l* A2 *QB1 [ [2] ] ] . 

MatrixExp [±* A3* QBl [ [3] ] ] .MatrixExp[l* AS *QB1 [ [4] ] ] . MatrixExp [1 * A3 * QBl [ [5] ]] . 
MatrixExp[i*A2*QBl[ [6]]] .MatrixExp [i * A3 * QBl [ [7] ] ] . MatrixExp [i * A8 * QBl [ [8] ] ] 

In[l4]:= UB2 : = MatrixExp[l * A3 * QB2 [ [1] ] ] .MatrixExp[l * A2 * QB2 [ [2] ] ] . 

MatrixExp [±* A3 *QB2[ [3] ] ] .MatrixExp[l* AS *QB2 [ [4] ] ] . MatrixExp [1 * A3 * QB2 [ [5] ] ] . 
MatrixExp[± * A2 * QB2 [ [6] ] ] .MatrixExp[l * A3 * QB2 [ [7] ] ] . MatrixExp [1 * A8 * QB2 [ [8] ] ] 
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In[15]:= (* The standard basis is chosen and the basis transformations are defined *) 



ln[16]:= 

ln[17]:= 
ln[18]:= 
ln[19]:= 
ln[20]:= 
ln[21]:= 



ln[25]:= 



I'l 












■('1 







.0. 






.1, 



eigal[k_] 
eigA2[k_] 
eigBl[k_] 
eigB2 [k_] 



Oal.v[ [k + 1] ] 

Oa2.v[ [k + 1] ] 

DBl.v[ [k + 1] ] 

DB2.v[ [k + 1] ] 



'•= eigAl[k] .ConjugateTranspose[eigAl [k] ] 
oA2 [k_] : = eigA2 [k] .ConjugateTranspose[eigA2 [k] ] 
oBl[l_] := eigBl[l] .ConjugateTranspose[eigBl [1] ] 
oB2[l_] := eigB2[l] .ConjugateTranspose[eigB2 [1] ] 

(* The Bell inequality 13 is defined *) 



ln[26]:= 



p[rho_, q_, w_, k_] := ^Tr [rho.KroneckerFroduct[q[Mod[ j + k, 3]], w[j]]] 



ln[27]:= i3[state_] := (p[state, oAl, oBl, 0] +p[state, oft2, oBl, -1] + 

p[state, oA2, oB2, 0] +p[state, oAl, oB2, 0]) - {p[state, oAl, oBl, -1] + 
p[state, oA2, oBl, 0] +p[state, oA2, oB2, -1] +p[state, oAl, oB2, 1]) 

ln[28]:= (* The maximization procedure is defined and the parameters are chosen *) 

ln[29]:= i3max[rho_] : = NMaximize[Re [i3 [rho] ] , domain, Maxlterations -» 200, 

Method -» { "NelderMead" , "ShrinkRatio" -» 0.8, "ContractRatio" -» 0.8, 

"Ref lectRatio" -» 1 . 6, "ExpandRatio" -» 1 . 6, "Tolerance" -» . 1, "RandomSeed" -» n} ] ; 

ln[30]:= (* The comnand 'compute' executes the 

maximization algorithm ten times with random simplices *) 

ln[31]:= compute [rho_] :=Do[max[n] = i3max[rho] , {n, 10}]; 

ln[32]:= (* The comnand 'maximum' yields the global 

maximum of 13 and the corresponding Euler angles *) 



ln[33]:= maximum := Array [max, 10] [ [Ordering[Array [max, 10], -1]]] 
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(* Example: QOO Bell state *) 



a = 2KroneckerProduct[v[ [ j + 1] ] , v[ [ j + 1] ] ] ; 

rhotest = S.ConjugateTranspose[S] 



/ 1 

3 











1 

3 











1 ' 

3 



















































































1 

3 











1 

3 











1 

3 



















































































1 

V 3 











1 

3 











1 

3 . 



compute [rhotest] 



maximum 

{{2.87293, {aAl[l] -> 5.72532, aAl[2] -> -0.520624, aAl[3] -> 1.21827, 



OAI [4] 


-> 


1 


71354, aAl[5] 


-> 


6.73672, 


cffil [6] 


-> 


9.07066, aAl[7] 


-> 1.79378, 


ctAl [8] 


-> 


1 


95697, ctA2[l] 


-> 


2.95736, 


ctA2 [2] 


-> 


3.32492, aA2[3] 


-> 0.0909413, 


aA2 [4] 


-> 


2 


28172, ca2[5] 


-> 


4.86642, 


ca2 [6] 


-> 


0.375456, aA2[7] 


-> 5. 67429, 


aA2 [8] 


-> 


3 


36199, aBl[l] 


-> 


5.06768, 


aBl [2] 


-> 


3.98113, aBl[3] 


-> 2.08775, 


aBl [4] 


-> 


1 


74692, aBl[5] 


-> 


3.09829, 


aBl [6] 


-> 


5.75545, aBl[7] 


-> 5.81317, 


oBl [ 8 ] 


-> 





711775, C(B2[1] 


-> 


2.40796 


, aB2[2] 


-> 


3.38539, aB2[3] 


-> 0.706663, 


aB2 [4] 


-> 


5 


1078, aB2[5] - 


> 4 


.34756, 


aB2[6] - 


> 2 


.80554, aB2[7] - 


> 4.09544, 


aB2 [ 8 ] 


-> 


9 


06912} } } 
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C MATHEMATICA notebook: Optimisation of Bj^ 



ln[1]:= (* The Gell-Mann matrices are defined *) 



ln[2]:= A2 ; 



( 
I 





A3 = 



A8 = 1/Sqrt[3] 



1 


10 








A5 : 



( 

I 




AlO : 



A15 = 1/Sqrt[6] 







-3 I 



ln[3]:= (* The generalized Euler angles for the four observables A1,A2,B1,B2 are introduced *) 
ln[4]:= QAl = Array [oAl, 15]; Element [QAl, Reals] 



Out[4]= (aAl(l) I aAl(2) | aAl(3) | aAl(4) | aAl(5) | aAl(6) | 

aAl(7) I q;A1(8) | q;A1(9) | q;A1(10) | q;A1(11) | q;A1(12) | q;A1(13) | q;A1(14) | q;A1(15)) ( 

ln[5]:= QA2 = Array [aA2, 15]; Element [QA2, Reals] 



Out[5]= (a;A2(l) | q;A2(2) | aA2(3) | aA2(4) | q;A2(5) | q;A2(6) | 

aA2(7) I q;A2(8) | q;A2(9) | q;A2(10) | q;A2(11) | q;A2(12) | q;A2(13) | q;A2(14) | q;A2(15)) 6 R 

ln[6]:= QBl = Array [aBl, 15]; Element [QBl, Reals] 

Out[6]= (a;Bl(l) I aBl(2) | aBl(3) | aBl(4) | aBl(5) | aBl(6) | 

aBl(7) I q;B1(8) | q;B1(9) I q;B1(10) | aBl(ll) | q;B1(12) | q:B1(13) | q;B1(14) | q;B1(15)) e R 

ln[7]:= QB2 = Array [aB2, 15]; Element [QB2, Reals] 

Out[7]= (a;B2(l) | q:B2(2) | q:B2(3) | q:B2(4) | aB2(5) | q;B2(6) | 

aB2(7) I q;B2(8) | q;B2(9) | q;B2(10) | aB2(ll) | q:B2(12) | q:B2(13) | q;B2(14) | q;B2(15)) e R 

ln[8]:= AV = Flatten [{QAl, QA2, QBl, £3B2}]; 

ln[9]:= Do[NB[i] = {AV[ [i] ] , 0, 2*Pi}, {i, 1, 60}]; domain = Array [NB, 60]; 
In[l0]:= (* The Euler angle parametrizations of UAl , UA2 , UBl , UB2 are defined *) 

ln[11]:= HAl := MatrixExp[i * A3 * QAl [ [1] ] ] . MatrixExp [ i * A2 * QAl [ [2] ] ] .MatrixExp[i * A3 * QAl [ [3] ] ] . 
MatrixExp[i* A5*QA1[ [4] ] ] .MatrixExp [i * A3 * QAl [ [5] ] ] . MatrixExp [i * AlO *nAl[ [6] ] ] . 
MatrixExp [i* A3* QAl [ [7] ] ] .MatrixExp [i * A2 * QAl [ [8] ] ] . MatrixExp [i * A3 * QAl [ [9] ] ] . 
MatrixExp [±*A5* QAl [[10]]] .MatrixExp[l * A3 * QAl [ [11] ] ] . MatrixExp [1 * A2 * QAl [ [12] ] ] . 
MatrixExp [±* A3* QAl [ [13]] ] . MatrixExp [1 * A8 * QAl [ [14] ] ] . MatrixExp [1 * Al 5 * QAl [ [15] ] ] 

ln[12]:= nA2 : = MatrixExp [1 * A3 * QA2 [ [1] ] ] .Matr±xExp[l* A2 *QA2 [ [2] ] ] . MatrixExp [1 * A3 * QA2 [ [3] ] ] . 
MatrixExp[±*A5*QA2[ [4]]] .MatrixExp[l* A3 *QA2 [ [5] ] ] . MatrixExp [1 * AlO *QA2 [ [6] ] ] . 
MatrixExp[i*A3*QA2[ [7] ] ] .MatrixExp[l* A2 *QA2 [ [8] ] ] . MatrixExp [1 * A3 * QA2 [ [9] ] ] . 
MatrixExp [i*A5*QA2[ [10] ] ] . MatrixExp [i * A3 * QA2 [ [11] ] ] . MatrixExp [i * A2 * QA2 [ [12] ] ] . 
MatrixExp [i * A3 *QA2 [ [13] ] ] . MatrixExp [i * A8 * QA2 [ [14] ] ] . MatrixExp [i * A15 *QA2 [ [15] ] ] 

ln[13]:= UBl := MatrixExp [1 * A3 * QBl [ [1] ] ] . MatrixExp [ 1 * A2 * QBl [ [2] ] ] . MatrixExp [1 * A3 * QBl [ [3] ] ] . 
Mat rixExp[±*A5*QBl[[4]]]. Mat r ixExp [1*A3*QB1[[5]]]. MatrixExp [1*A10*QB1[[6]]]. 
MatrixExp[±*A3*QBl[ [7]]] .MatrixExp [1 * A2 * QBl [ [8] ] ] . MatrixExp [1 * A3 * QBl [ [9] ] ] . 
MatrixExp[i*A5*QBl[ [lOj] ] .MatrixExp[l* A3 *QB1 [ [11] ] ] . MatrixExp [1 * A2 * QBl [ [12] ] ] . 
MatrixExp [±* A3* QBl [ [13] ] ] . MatrixExp [1 * A8 * QBl [ [14] ] ] . MatrixExp [1 * Al 5 * QBl [ [15] ] ] 

ln[14]:= HB2 := MatrixExp [i * A3 *QB2 [ [1] ] ] . MatrixExp [ i * A2 * QB2 [ [2] ] ] . MatrixExp [i * A3 * QB2 [ [3] ] ] . 
MatrixExp[±*A5*QB2[ [4]]] .MatrixExp [1 * A3 * QB2 [ [5] ] ] . MatrixExp [1 * AlO *QB2 [ [6] ] ] . 
MatrixExp [±* A3 *QB2[ [7] ] ] .MatrixExp[l* A2 *QB2 [ [8] ] ] . MatrixExp [1 * A3 * QB2 [ [9] ] ] . 
MatrixExp [±*A5*QB2[ [10] ] ] .MatrixExp[l* A3 *QB2 [ [11] ] ] . MatrixExp [1 * A2 * QB2 [ [12] ] ] . 
MatrixExp [±* A3* QB2[ [13] ] ] . MatrixExp [1 * A8 * QB2 [ [14] ] ] . MatrixExp [1 * Al 5 *QB2 [ [15] ] ] 
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In[l5]:= (* The standard basis is chosen and the basis transformations are defined *) 



ln[16]:= 

ln[17]:= 
ln[18]:= 
ln[19]:= 
ln[20]:= 
ln[21]:= 

ln[25]:= 



1 








'"1 




(0\ 







1 






















1 







.0, 




,0 




,0 




1 



eigAl[k_] 
eigA2[k_] 
eigBl[k_] 
eigB2 [k_] 



= UAl.v[ [k + 1] ] 

= UA2.v[ [k + 1] ] 

= UBl.v[ [k + 1] ] 

= UB2.v[ [k + 1] ] 



oAl [k_] : = eigAl [k] . Con jugateTranspose [eigAl [k] ] 

oA2 [k_] : = eigA2 [k] .Con jugateTranspose [eigA2 [k] ] 
oBl : = eigBl [1] . Con jugateTranspose [eigBl [1] ] 

oB2 [1_] : = eigB2 [1] . Con jugateTranspose [eigB2 [1] ] 

(* The Bell inequality 14 is defined *) 



ln[26]:= p[rho_, q_, w_, k_] := ^Tr [rho.KroneckerProduct[q[Mod[ j + k, 4]], w[j]]] 
ln[27]:= i4[state_] : = 

Sum[ (1 - 2 *k/ 3) ((p[state, oAl, oBl, k] +p[state, oA2, oBl, -1-k] +p[state, oA2, oB2, k] + 
p[state, oAl, oB2, -k]) - (p[state, oAl, oBl, -1-k] +p[state, oA2, oBl, k] + 
p[state, oA2, oB2, -1-k] +p[state, oAl, oB2, 1 + k] ) ) , (k, 0, 1}] 

ln[28]:= (* The maximization procedure is defined and the parameters are chosen *) 

ln[29]:= i4inax[rho_] : = NMaximize[Re [i4 [rho] ] , domain, Maxlterations -> 200, 

Methods { "NelderMead" , "ShrinkRatio" ^0.8, "ContractRatio" ^0.8, 

"Ref lectRatio" -» 1 . 6, "ExpandRatio" -» 1 . 6, "Tolerance" -» . 1, "RandomSeed" -» n} ] ; 

ln[30]:= (* The command 'compute' executes the 

maximization procedure ten times with random simplices *) 

ln[31]:= confute [rho_] :=Do[max[n] = i4max[rho] , {n, 10}]; 

ln[32]:= (* The conmand 'maximum' yields the global 

maximum of 14 and the corresponding Euler angles *) 



ln[33]:= maximum := Array [max, 10] [ [Ordering[Array [max, 10], -1]]] 
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(* Example: QOO Bell state *) 



S = - ^KroneckerProduct[v[ [ j + 1] ] , v[ [ j + 1] ] ] ; 

^ 3=0 



rhotest = S.ConjugateTranspose[S] 
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compute [rhotest] 
maximum 

{{2. 89624, {ctAl[l] -> 3.07265, ctAl[2] -> 6.15285, ctAl [3] -> 5. 67497, ctAl[4] -> 1. 92491, 
cffil[5] -> 2 .05345, C(A1[6] -> 3.86138, ctAl[7] -> 2.57146, C(A1[8] -> 1 . 96996, 
ciAl[9] -> 5.06114, aAl[10] -> 2.7175, aAl[ll] -> 2.73205, aAl[12] -> 4.04356, 
CIA1[13] -> 3.52083, C(A1[14] -> 3.34487, C(A1[15] -> 2.0385, C(A2[1] -> 0.72059, 
C(A2[2] -> 3.20978, ca2 [3] -> 3.40004, aA2[4] -> 2.25048, ca2 [5] -> 4.36813, 
aA2[6] -> 3.47318, aA2[7] -> 5.05212, aA2r8] -> 6.39357, aA2[9] -> 2.44141, 
CIA2[10] -> 2.28185, aA2[ll] -> 3.52065, aA2:i2] -> 3.93748, aA2[13] -> 1.8451, 
aA2[14] -> 1.67687, C(A2[15] -> 4.53126, aBl[l] -> 5.23488, aBl[2] -> 1.76584, 
aBl[3] -> 2 .8723, aBl[4] -> 4. 64033, C(B1[5] -> 4.4815, C(B1[6] -> 2.10096, 
aBl[7] -> 2 .32368, aBl[8] -> 4.54152, c(Bl[9] -> 3. 61002, aBl[10] -> 1.1092, 
aBl[ll] -> 4.20415, aBl[12] -> 1.68768, aBl[13] -> 2.77742, aBl[14] -> 2.38235, 
aBl[15] -> 2.27769, C(B2[1] -> 1.55276, aB2[2] -> -0.00921043, aB2[3] -> 2 .56765, 
C(B2[4] -> 5.51796, C(B2 [5] -> 3. 68283, C(B2[6] -> 2.72437, C(B2[7] -> 4.73605, 
aB2[8] -> 4 . 09868, aB2[9] -> 2.13773, aB2[10] -> 2.89167, C(B2[11] -> 1.43804, 
aB2[12] -> 3. 99496, C(B2[13] -> 1.30861, C(B2[14] -> 3. 90506, C(B2[15] -> 2 . 92697 })} 
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D Computation details on p = ^—^—^1 + aPo^ + pPifi 



PPT: 



w — e 
co,o = 

Cl,0 = 
Ck,l = 

' detBg — 



2m/3 



1 


— q; 


-P 




9 




1 


— a 


-P 




9 




1 


— a 




9 



+ a 



:,k{S-t) 



V(A;,0^{(0,0),(1,0)} 



^ Ck,s+tW 
k=0 

( i(2Q; + 2/3 + 1) 




729 



\(a + ^e^^ |(-q;-/3 + 1) ^ 
(2a + 2/3 + l)(8a^ + 8/3^ - lla/3 + 2a + 2/3 - 1) 
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E Computation details on p = ^ ^ 1 + aPo,o + (^Pi,o + 7^2,0 



PPT: 



w — e 



2m/3 



co,o 

Cl,0 
C2,0 
Ck,l 

► detBn — 



1 


— a 




-7 






9 




1 


— a 




-7 






9 




1 


— a 


-P 


-7 






9 




1 


— a 


-/5 


-7 



+ q; 
+ P 

+ 7 



Ck,s+tW 



k{s-t) 



k=0 



V(fc,0^{(0,0),(l,0),(2,0)} 



/ i(2a + 2/5 + 27 + 1) 



9 




1 



i(-a-/9-7 + l) I (a + pe-^ + -fe^ 
l(a + pe^ + -fe-^^ i(-a-/?-7 + l) 



(2q; + 2/3 + 27 + 1) 



729 

{8a^ + 8/52 + 87^ + 2a + 2/5 + 27 - ll/3a - 11^7 - II/37 - 1) 
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F Computation details on p = ^-^-^-^1 + aPo,o + (^Pi,q + 7^o,i 



PPT: 



w — e 



co,o 

Cl,0 

Co,i 

^^0 



27ri/3 



1 


— a 


-(3- 


7 






9 




1 


— a 


-P- 


7 






9 




1 


— a 


-(5- 


7 






9 




1 


— a 




7 



+ q; 
+ /? 

+ 7 



if 



k{s-t) 



V (A;, 0^ {(0,0), (1,0), (0,1)} 



fc=0 

/ |(2a + 2/5-7 + l) I 

2 i(-a-/3-7 + l) i(^a + /3e- 







\(a + l3e^^ i(-a-/?+27 + l) ^ 



deiSo = -^(-16q;^ - 16/3^ - 167=^ + 6/3q;^ + 670;^ + 67^0; + 6/3^0; + 6/3^7 
729 

,,2 1 o„,2 1 o /32 1 0^,2 



+ 6/37^ - 12a^ - 12/3^ - 127^ + 3/3a + 870; + 3/37 - 15/37a + 1) 
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